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The  aerodynamic  characteristics  of  thirteen  distinct  chaff  dipole  con- 
figurations were  determined  from  drop  tests  of  individual  elements  in  a 
special  enclosed  test  chamber.  The  dipole  motion  and  trajectory  were 
recorded  by  multi-image  photographs  taken  by  orthogonal  still  cameras 
equipped  with  specially  designed  synchronized  rotating  shutters. 

The  dynamic  behavior  and  descent  rate  of  the  dipoles  was  found  to  (cont'd 
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depend  greatly  upon  the  principal  cross-section  dimensions  of  the  filaments. 

Aerodynamic  forces  and  moment  coefficients  for  each  dipole  configuration 
were  computed  from  the  photographic  multi-image  motion  data  using  photo- 
grammetric  and  aerodynamic  data  reduction  programs  developed  as  part  of 
the  effort.  The  resulting  aerodynamic  coefficient  data  were  successfully  cor- 
related with  Reynolds  number,  angle  of  attack,  and  various  other  parameters. 
While  the  force  coefficients  were  found  to  be  large  and  in  general  agreement 
with  the  various  theories  for  creeping  flow,  the  moment  coefficients  were 
extremely  small  and  resulted  primarily  from  configurational  asymmetries. 

Using  both  the  experimental  data  and  theory,  aerodynamic  coefficient 
tables  for  representative  dipole  configurations,  suitable  for  6-DOF  simulation 
of  dipole  motion,  were  prepared.  These  aerodynamic  data  were  subsequently 
used  in  conjunction  with  a 6-DOF  Monte  Carlo  trajectory  program,  modified 
for  inclusion  of  stochastic  atmospheric  turbulence,  for  preliminary  simulation 
of  chaff  dipole  motion  in  both  quiescent  and  turbulent  atmospheres. 

Turbulence  was  found  to  have  a large  effect  on  the  translational  motion  of 
the  dipole,  but  only  a small  effect  on  its  angular  motion. 
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SECTION  I 


INTRODUCTION 


Chaff,  as  a broad  band  passive  electromagnetic  countermeasure, 
has  found  extensive  use  since  World  War  II.  Apart  from  slip  coating  to 
help  separation  and  strip  coating  to  assure  different  orientations  under 
free-fall,  no  systematic  effort  has  beer,  made  to  improve  the  perform- 
ance of  atmospheric  chaff  clouds  taking  into  account  the  unique  flight 
performance  characteristics  of  the  elements  themselves.  Before  chaff 
improvements  can  be  initiated  in  a logical  and  orderly  manner,  the 
dynamic  behavior  of  the  chaff  dipoles  must  be  fully  understood  and  the 
means  for  predicting  the  dipole  trajectory  and  angular  motion  during 
both  short  and  long  time  spans  must  also  be  available. 

Chaff  dipoles,  because  of  their  extremely  small  cross-sectional 
dimensions  and  great  slenderness,  possess  unusual  flight  characteristics 
The  flight  Reynolds  number,  which  is  of  the  order  of  unity,  approaches 
the  Stokes  flow  regime  at  the  lowest  Reynolds  number,  and  the  transi- 
tional regime  where  fluid  vortices  are  formed  at  the  largest  flight 
Reynolds  numbers.  The  extremely  slow  descent  rate  of  chaff  elements 
(which  is  of  the  order  of  1 FPS)  makes  them  highly  sensitive  to  small 
air  currents  and  turbulence. 

The  present  effort  was  undertaken  for  the  purpose  of  I)  detailed 
measurement  of  the  motion  and  aerodynamic  coefficients  of  individual 
chaff  dipoles  provided  by  the  Air  Force  Avionics  Laboratory, 

2)  determining  the  effect  of  dipole  configuration  on  the  characteristic 
aerodynamic  parameters,  and  3)  for  development  of  a 6-DOF  (six- 
degrees-of-freedom)  trajectory  simulation  model  for  prediction  of 
dipole  motion  in  the  presence  of  realistic  winds  and  turbulence. 

A prior  study  of  chaff  cloud  dynamics,  as  affected  by  the  dipole 
aerodynamic  characteristics,  is  described  in  Reference  1.  Although 
Reference  1 reports  on  the  aerodynamic  properties  of  a single  repre- 
sentative foil-type  dipole,  insufficient  aerodynamic  coefficient  data  were 
determined  to  permit  6-DOF  type  motion  simulations. 


Avrin,  J.  S.  it  M.  V.  McCafferty,  The  KMS  Technology  Center,  "Chaff 
Cloud  Dynamics,"  Air  Force  Avionics  Laboratory,  Wright-Patterson 
AFB,  Ohio,  Report  No.  AFAL-TR-73-286,  September  1973. 
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The  experimental  measurement  of  dipole  aerodynamic  character- 
istics presents  several  problems  both  as  a result  of  the  small  dipole  size 
(minimum  diameter  of  0.001  inches)  and  the  sensitivity  of  the  dipole 
motion  to  small  air  currents.  Because  of  t ese  factors  it  was  concluded 
at  the  outset  of  the  program  that  special  test  facilities,  test  procedures, 
and  photographic  techniques  would  be  required  to  accurately  measure  the 
chaff  motion  and  the  associated  aerodynamic  fore ec  and  moments.  As  a 
result,  a unique  fully-enclosed  drop  test  chamber  was  fabricated,  with 
an  integral  lighting  and  photographic  documentation  system. 

To  document  the  dipole  motion,  three  orthogonally  aligned  still 
cameras,  each  equipped  with  a specially  designed  synchronous  rotating 
shutter,  was  provided  in  three  cells  attached  to  the  main  drop  chamber. 
This  method  of  multi-image  photography,  using  single  film  negatives, 
was  preferred  over  conventional  motion  pictures  because  of  the  greater 
flexibility,  simplified  data  processing,  and  improved  accuracy. 

Satisfactory  illumination  of  the  dipoles  in  flight  was  a significant 
problem  and  required  considerable  development  effort  with  regard  to 
lighting,  baffling,  exposure  control,  and  film  selection  and  processing. 

A total  of  262  separate  documentary  tests  were  accomplished  in 
the  aerodynamic  test  chamber,  using  13  dipole  configurations  differing 
in  construction,  cross-section,  and  length. 

Aerodynamic  data  reduction  and  analysis  were  accomplished  in 
two  phases.  In  the  first  phase  the  dipole  trajectory  and  attitude  was 
determined  as  a function  of  time.  In  the  second  phase  the  motion 
parameters,  including  all  velocity  and  acceleration  components,  were 
obtained  by  moving  polynomial-arc  smoothing  and  differentiation  formu- 
las, and  subsequently  the  total  aerodynamic  force  and  moment  coeffi- 
cients were  computed  from  the  solution  of  the  applicable  equations  of 
motion.  The  aerodynamic  coefficient  dependence  upon  angle-of-attack 
and  Reynolds  number  was  determined  by  data  correlation  techniques. 

To  complement  the  experimental  aerodynamic  analysis  a 
comprehensive  theoretical  analysis  of  chaff  dipole  aerodynamics  was 
undertaken.  These  analyses  provided  a means  for  estimating  the  effects 
of  dipole  roll  orientation,  configurational  asymmetry  and  aerodynamic 
damping,  which  could  not  be  determined  experimentally. 

A significant  part  of  the  present  effort  was  devoted  to  the  develop- 
ment of  a 6- DOF  computer  program  and  simulation  capability.  This  was 
accomplished  by  adaptation  of  an  existing  comprehensive  6- DOF  program 
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Preliminary  6 -DOF  motion  simulations  were  accomplished  for 
r.pr««  J“po,,  configurations  to  iilustrate  dipoic  tnotton  char- 
acteristics  in  both  quiescent  and  turbulent  atmospheres. 
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SECTION  II 


TEST  FACILITIES 


I.  Enclosed  Test  Chamber 


All  instrumented  drop  tests  of  chaff  dipoles  were  accomplished  in 
a special  aerodynamic  test  chamber  as  depicted  in  Figure  1.  The  test 
chamber  provides  an  isolated  atmospheric  environment  (at  approximately 
sea  level  atmospheric  pressure),  free  of  air  currents  in  the  surrounding 
laboratory,  and  is  constructed  in  such  a manner  as  to  enhance  the  illum- 
ination and  photographic  documentation  of  single  dipoles  in  free  descent. 

The  primary  viewing  chamber  is  approximately  a 4-ft  cube,  which 
is  fully  illuminated  by  twelve  500-W  photoflood  bulbs  in  12-inch  alumin- 
ized reflectors  placed  around  the  periphery  as  indicated  by  Figure  1. 
Contiguous  to  the  primary  chamber  are  three  cells  for  mounting  of  three 
fixed  motion-recording  cameras  and  two  auxiliary  chambers,  which 
extend  the  background  viewed  by  the  two  horizontally  aligned  cameras. 

The  background  extensions  prevent  light  from  being  reflected  off  the 
chamber  surfaces  and  onto  the  camera  lenses,  such  that  a high  contrast 
is  maintained  between  the  illuminated  dipoles  and  the  background.  To 
further  improve  the  contrast  all  of  the  interior  surfaces  of  the  test 
chamber  are  painted  flat  black  and  the  background  walls  in  view  of  the 
cameras  are  draped  with  black  velour.  In  addition,  a series  of  primary 
and  secondary  baffles  are  provided  to  further  control  the  background 
illumination. 

The  cell  directly  above  the  test  chamber  contains  the  chaff  release 
mechanisms  as  well  as  the  downward  viewing  camera.  The  overall  height 
of  the  chamber  is  10  ft.  , thus  permitting  dipole  free-fall  of  up  to  six  feet 
prior  to  entry  of  the  viewing  section.  The  overall  length  and  width  of  the 
test  chamber  is  16  ft.  by  16  ft. 

The  size  of  the  primary  viewing  chamber  was  dictated  by  three 
considerations.  First,  the  internal  volume  and  dimensions  had  to  be 
adequate  to  eliminate  wall  effects  on  the  chaff  aerodynamics.  Secondly, 
the  chamber  air  mass  had  to  be  sufficiently  Bmall  that  it  would  damp  out 
air  movements  rapidly.  Third,  the  chamber  size  had  to  be  small  enough 
to  allow  a very  high  level  of  dipole  illumination,  since  light  intensity 
decreases  with  distance  squared. 

The  10-ft.  maximum  drop  I'.eight  is  quite  ample  for  investigation  of 
dipole  steady- state  descent.  The  approximate  distance  to  reach  terminal 
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Figure  1.  Chaff  Aerodynamic  Test  Chamber 
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descent  velocity  when  a dipole  is  released  with  no  initial  velocity  is 
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V88  = steady- state  velocity 
t = time 


a = 0,99  (i.  e. , time  to  achieve  0.99  V88). 


Using  the  fastest  terminal  descent  velocity  observed  in  preliminary  tests 
(Vss  = 3 ft/sec).  Equation  (1)  results  in  a transient  descent  distance  of 
0.55  ft.  , which  occurs  in  only  0.25  second. 


2.  Instrumentation  and  Equipment 

Cameras  Three  Calumet  * in.  x 5 in.  view  cameras  with 
210  mm  f/6.3  Caltar  lenses  are  used  to  record  the  dipole  motion.  The 
cameras  are  fixed  in  position  and  alignment.  Each  camera  station  is 
provided  with  an  independently- supported  special  rotating- shutter 
system  to  "strobe"  the  dipole  images.  The  integral  shutter  in  each  lens 
assembly  is  also  retained  for  initiating  and  terminating  the  motion 
sequences.  Both  shutters  are  electrically  powered,  synchronized,  and 
remotely  operable.  Figure  2 depicts  a camera  and  rotating  shutter 
installation. 

The  cameras  are  positioned  and  aligned  such  that  the  lens  axes  of 
cameras  1 and  2 are  both  horizontal  and  orthogonal.  Camera  - 3 is 
positioned  above  the  viewing  chamber  and  the  lens  axis  is  inclined  16.6 
degrees  from  vertical,  such  that  the  dipole  images  during  vertical 
descent  are  separated.  Cameras  l and  3 constitute  a second  orthogonal 
pair. 


Rotating  Shutter  System  Each  camera  station  is  equipped  with 
a secondary  rotating-disk  shutter  providing  intermittent  strobing  at  a rate 
of  60  images/sec  (or  optionally,  120  images/sec  with  two  slot  openings) 
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Figure  2.  Camera  System  and  Rotating  Shutter 


and  an  exposure  duration  of  1/1000  second.  These  shutters  were 
specially  designed  for  the  chaff  aerodynamic  test  program. 

Each  shutter  is  comprised  of  a 12-inch  diameter,  0.030-inch 
thick  aluminum  disk,  driven  by  a 1/50  horsepower  3600  rpm  hysteresis 
synchronous  motor.  An  AC  generator  is  attached  to  each  motor  to 
determine  mechanical  phase  orientation  and  to  facilitate  synchronization. 
The  disk  and  motor  assembly  are  installed  in  a close  fitting  housing, 
which  acts  ar  a light  baffle  and  as  a means  for  reducing  both  the  aero- 
dynamic frictional  torque  on  the  disk  and  induced  air  currents  within  the 
test  chamber. 

Fig  ire  3 shows  a schematic  of  the  rotating  shutter  system.  The 
three  rotating  shutters  synchronized  at  the  beginning  of  each  testing 
period  by  mechanically  varying  the  disk  phase  position.  The  phase  can  be 
maintained  for  long  periods  of  time  and  in  normal  operation  the  shutters 
are  aligned  each  morning  and  left  running  continually  throughout  the  day. 

Chaff  Release  Mechanism  The  device  used  to  hold  and  release 
the  individual  chaff  dipoies  is  depicted  in  Figure  4.  The  holder  is  re- 
motely operable  and  can  be  pre-positionea  to  the  desired  dipole  release 
orientation. 

Pairs  of  0.008-inch  stainless-steel  wires  with  a tweezer  action 
holds  the  dipole  until  a solenoid  actuated  rod  displaces  the  dipole  from 
the  holder.  The  solenoid  has  a nominal  stroke  of  0.  1 inch  with  0.1  oz. 
of  force  at  this  displacement.  An  initial  velocity  (slightly  less  than  the 
terminal  velocity)  is  usually  imparted  by  the  release  mechanism;  how- 
ever, the  release  velocity  can  be  varied  by  changing  the  applied  voltage 
to  the  solenoid. 

A bolder  assembly  with  a cluster  of  six  identical  release  mechan- 
isms is  shown  in  the  photographs,  Figure  5.  The  complete  holder 
assembly  and  wiring  harness  are  attached  to  a strut  which  spans  the 
test  chamber.  This  strut  can  be  positioned  vertically  at  any  desired 
height  within  the  chamber. 

3.  MCAS  Drop  Facility 

Test  drops  of  chaff  clusters  (to  determine  dip  rsive  character- 
istics at  drop  heights  up  to  150  ft.  ) are  accomplish  cd  in  the  Marine 
Corps  Air  Station  (MCAS)  airship  hangar  at  Santa  Ana,  Calif.  This 
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Figure  4.  Filament  Holder /Release  Mechanism 


hangar  has  a height  of  155  ft.  and  floor  space  of  300  ft.  x 1100  ft.  The 
facility  provides  a relatively  quiescent  atmospheric  environment  during 
the  night  when  ail  of  the  hangar  doors  and  entrances  are  closed. 

Chaff  cluster  drops  are  accomplished  by  use  of  a trap*door  type 
holder  suspended  from  the  hangar  ceiling  and  positioned  by  a variable 
length  tether.  The  surface  of  the  trap  door  is  teflon  coated  to  prevent 
the  dipoles  from  sticking. 
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SECTION  III 


TEST  ARTICLES 


Test  dipoles  were  selected  from  four  basic  types  of  chaff 
supplied  by  the  Air  Force  project  monitor.  Table  1 summarizes  the 
chaff  type  specifications  and  designates  the  dipole  sizes  selected 
for  the  aerodynamic  test  program.  Articles  1 through  11  are  those 
originally  selected  for  the  test  program,  while  articles  12  and 
13  were  Included  for  the  purpose  of  extending  the  scope  of  the 
aerodynamic  investigations. 

Representative  dipoles  have  been  examined  in  great  detail. 

It  was  noted  that  the  aluminum  coated  glass  dipoles  exhibit  con- 
siderable surface  roughness.  Likewise,  the  edges  of  the  foil  chaff 
exhibit  email  irregularities.  The  nominal  width  (or  diameter)  of 
each  dipole  type  was  approximately  verified  by  using  a microscope 
in  conjunction  with  an  optical  comparator. 

Several  dipole  configurations  Incorporate  a longitudinal 
V-bend  for  increased  stiffness  (see  Table  1) . The  V-bend  shape 
varies  with  chaff  type,  and  is  qualitatively  described  by  the 
following  sketch. 


Type  RR- 141 /AL  Type  RR- 149 /AL  Type  RR-39A/AL 

V*Bend  Configurations 

The  nominal  chaff  dimensions  are  depicted  by  the  following  sketch 
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Table  1.  DESIGNATION  OF  CHAFF  TYPES 


"0  TO 
o>  u 


<2  >.  >. 
y TO  T3  TO 

l § 3 S 

« i-l  i-J  >-5 


tj 

TO 

TO 

C 

•0 

c 

T3 

TO 

TO 

c 

C 

0 

0) 

C 

e 

c 

0) 

0) 

A 

A 

V 

0) 

01 

A 

A 

> 

i 

XI 

XI 

i 

• 

> 

> 

i 

1 

1 

> 

> 

0 

o 

> 

> 

> 

c 

c 

a a 


TO  TO 

c c 

0)  V 

A A 

i i 

> > 


in 

eg 

oo 

00 

N 

A 

N 

in 

NO 

CO 

vO 

CO 

H 

in 

vO 

M 

« 

P-X 

• 

o 

• 

o 

• 

«-4 

• 

• 

o 

• 

• 

o 

X 

X 

X 

X 

X 

X 

X 

X 

X 

<NJ 

N 

r~ 

in 

in 

m 

in 

in 

in 

in 

in 

in 

X 

X 

X 

X 

co 

CO 

m 

CO 

•# 

n}< 

x* 

nf 

o 

o 

o 

o 

o 

O 

o 

o 

O 

•«* 

• H 

•w* 

o 

o 

o 

o 

o 

o 

o 

o 

O 

£ 

£ 

£ 

£ 

o 

o 

o 

o 

o 

o 

o 

o 

O 

eg 

^4 

•— 4 

X 

X 

X 

X 

X 

X 

X 

X 

X 

00 

00 

NO 

NO 

A 

o 

o 

o 

o 

o 

o 

o 

o 

© 

o 

• 

o 

• 

O 

• 

o 

• 

o 

• 

o 

• 

O 

• 

© 

• 

o 

• 

— i vO  00  <\] 

Xfc  =*:  -h 

_ =*: 


® — c 
**  60 
a >h 

< “ 
^ 01 

Q 


o n n « 


■ 1 


Note:  1)  All  items  slip  coated. 

2)  All  dimensions  in  inches. 


in  the  case  of  dipoles  with  V-bend  the  nominal  width  corresponds  to  the 
equivalent  flat  strip. 

Using  the  above  nominal  dimensions,  the  weight  and  inertia  of 
each  type  dipole  was  computed  based  on  the  following  material  densities. 


Material 

fiberglass 

aluminum  foil  or  coating 
slip  coating 


Density 
0.075  lb/in3 
0.0975  lb/in3 
0.05  oz/1000  in^ 


The  computed  dipole  weights  are  summarized  in  Table  2.  To 
verify  the  calculated  weights,  actual  average  weights  for  several  repre- 
sentative dipole  types  were  determined  with  a microbalance.  Although 
actual  weights  varied  from  sample  to  sample,  relatively  good  agreement 
with  the  calculated  weights  was  achieved  for  all  of  the  designated  chaff 
types  except  the  2 mil  x 2 inch  dipole,  which  was  found  to  have  an  actual 
average  weight  of  only  3.76  x 10' ' lbs,  about  60  percent  of  the  nominal 
calculated  value. 

Measured  weight  data  were  also  determined  for  the  RR-39A/AL 
specification  chaff,  because  these  dipoles  have  an  additional  lead  strip 
coating  which  is  used  to  vary  the  center-of-mass.  For  example,  the 
weight  of  the  .040  x .00045  x 1.12  dipole  varies  from  about  2.07  - 
2.39  x 10'^  lbs.  depending  upon  the  amount  of  atrip  coating.  Average 
weights  for  test  articles  13  and  14  are  2.22  x 10'^  and  0.960  x 10”^  lbs, 
respectively. 
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CHAFF  DIPOLE  WEIGHT  AND  INERTIA  SUMMARY 
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SECTION  IV 


TEST  CONDITIONS  AND  PROCEDURES 


The  test  program  utilizing  the  special  test  chamber  was  divided 
into  the  following  major  phases: 

1)  photogrammetric  calibration 

2)  development  of  photographic  procedure 

3)  investigation  of  test  chamber  turbulence  and  air  currents 

4)  chaff  dipole  steady-state  descent  tests 

5)  chaff  dipole  transient  motion  tests 

6)  special  investigations  of  chaff  dynamics. 

A total  of  423  tests  were  accomplished  of  which  161  were  for 
development  and  calibration  of  the  photogrammetric  system  and  262  for 
documentation  of  chaff  dipole  flight. 

1.  Photogrammetric  Calibration 

The  photogrammetric  calibration  procedures  are  described  in 
Appendix  A.  The  final  reference  system  for  measurement  of  dipole 
position  and  attitude  consisted  of  a single  plumb  wire  with  three  bead- 
type  reference  marks,  each  12  inches  apart.  The  plumb  line  and  beads 
are  clearly  visible  in  the  photographs  taken  from  each  of  the  three 
cameras,  but  are  so  located  that  the  dipole  motion  is  never  obscured. 

2.  Development  of  Photographic  Procedure 

Special  photographic  procedures  were  required  in  order  to 
obtain  distinguishable  chaff  images  on  film.  At  all  times  there  was  a 
compromise  between  high  lens  f stops  to  realize  depth  of  field  and 
lower  f stops  to  improve  exposure. 

The  exposure  time  of  0.001  sec.  (provided  by  the  slot  geometry 
of  rotating  shutter)  was  chosen  because  it  results  in  an  image  size  which 
can  be  resolved  on  film.  The  relationship  between  image  size  and  reso- 
lution is  shown  in  Figure  7.  The  various  sizes  shown  in  Figure  7 relate 
to  negative  film  resolution  in  such  a way  that  one  unit  corresponds  to  100 
line  parts  per  millimeter  (Lppmm).  The  nominal  filament  diameter  is 
shown  to  be  significantly  smaller  than  the  resolving  capability  of  the  film. 


but  with  the  streaking  provided  by  the  rotating  shutter  the  image  size  is 
significantly  larger  than  the  resolution.  Also,  the  24-inch  out-of-focus 
image  width  of  approximately  25  resolutions  is  readable  on  the  negatives 
and  prints  because  the  circle  of  confusion  is  0.01  inches.  Thus,  with 
adequate  illumination  a 48-inch  depth  of  field  (corresponding  to  the  test 
chamber  dimensions)  is  achievable. 


The  film  finally  selected  for  the  documentary  tests  was  Kodak 
Royal-X  pan  film,  Type  4166,  in  4 in.  x 5 in.  Estar  sheets.  It  has  a 
rated  ASA  index  of  4000,  which  is  the  highest  commercially  available 
speed  with  standard  processing.  Further  improvement  in  photographic 
image  quality  was  achieved  by  pre-exposing  each  negative  to  a density  of 
0.35  to  0.45  (4-percent  opaque)  and  by  extended  film  developmen  (15 
minutes)  with  Kodak  D-19  developer.  The  pre-exposure  was  necessary 
to  force  the  emulsion  above  its  threshold  value,  and  resulted  in  a- 
effective  speed  higher  by  several  f stops. 


The  photographic  system  was  also  required  to  produce  up  to  180 
dipole  images  per  film  negative  without  background  over -illumination. 
This  number  of  images  represents  three  seconds  of  flight  when  the 
rotating  shutter  is  operating  at  60  images/second.  The  directly  illum- 
inated dipoles  have  a probable  reflection  index  of  0.98,  therefore,  if  a 
factor  of  10  difference  in  background  reflectance  is  necessary  to  ade- 
quately perceive  an  image,  the  allowable  background  reflectance  is 
approximately 


0.  98 
10  x 180 


0.005 


This  low  level  of  reflection  was  achieved  by  use  of  background  materials 
having  a reflectance  of  less  than  0.02  and  by  locating  the  background 
cells  and  baffles  such  that  all  visible  background  light  rays  were  doubly 
reflected. 


Image  correlation  between  the  three  camera  negatives  was  accom- 
plished in  one  of  four  ways. 

1)  dipole  initial  separation  from  the  holder. 

2)  first  image  at  photoflood  turn  on 

3)  last  image  at  photoflood  turn  off 

4)  vertical  position  correlation. 
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3.  Investigation  of  Test  Chamber  Air  Currents 

The  large  quantity  of  energy  placed  into  the  test  cell  by  the  photo- 
flood lamps  was  found  to  produce  a noticeable  convective  movement  of 
air.  The  onset  of  this  convective  movement  was  investigated  by  injection 
of  a white  talc  cloud  prior  to  lamp  turn-on.  Observations  of  the  talc- 
cloud  showed  that  the  air  mass  in  the  main  4-ft.  cube  begins  a vertical 
rise  closest  to  the  lights,  commencing  about  5 seconds  after  turn-on, 
and  extends  into  the  center  of  the  test  chamber  after  10  to  15  seconds. 
Therefore,  in  the  aerodynamic  tests,  the  photof.ood  lamps  are  turned 
on  either  at  the  instant  of  chaff  release  or  at  the  approximate  time  of 
chaff  entry  into  the  viewing  chamber. 

The  atmospheric  temperature  rise  during  each  test  was  esti- 
mated to  be  about  0.06“  F,  due  primarily  to  water  vapor  absorption. 

In  addition,  a significant  increase  in  wall  temperature  was  observed. 

To  maintain  thermal  equilibrium  and  to  provide  for  damping  of  all  con- 
vective air  currents,  a minimum  of  10  minutes  was  allowed  between 
all  drop  tests. 

4.  Dipole  Steady-State  Descent  Tests 

Aerodynamic  tests  to  determine  the  dipole  steady- state  descent 
characteristics  were  accomplished  with  the  release-holder  mechanism 
positioned  from  two  to  four  feet  above  the  photographic  viewing  chamber, 
thus  permitting  the  dipoles  to  attain  terminal  velocity  prior  to  entering 
the  viewing  chamber.  One  dipole  (or  in  some  instances  two  dipoles  of 
widely  different  length)  were  released  and  followed  visually  until  they 
reached  the  top  of  the  photographic  viewing  chamber  and  at  this  time  the 
photoflood  lights  were  activated.  All  steady-state  tests  were  accom- 
plished with  the  release  holder  mechanism  in  the  horizontal  orientation. 

Test  Numbers  162  - 270  and  388  - 390  were  for  the  purpose  of 
documenting  chaff  steady- state  flight  behavior.  These  tests  were  per- 
formed using  all  of  the  test  articlee  designated  in  Table  2 (test  articles 
1 through  13).  Each  dipole  type  was  dropped  several  times  to  insure 
adequate  photographic  coverage. 

5.  Dipole  Transient  Motion  Tests 


Tests  to  determine  the  dipole  transient  motion  were  accomplished 
with  the  release-holder  mechanism  positioned  at  the  top  of  the  4-ft.  cube 
test  chamber,  in  view  of  all  three  cameras.  The  chaff  dipole  motion  was 


recorded  from  the  instant  of  release.  Four  different  orientations  of  the 
release  mechanism  were  used  to  obtain  initial  chaff  dipole  attitudes  of 

1)  horizontal 

2)  45-degrees  from  horizontal 

3)  60-degrees  from  horizontal 

4)  vertical. 

Test  Numbers  271  - 387  were  for  the  purpose  of  recording  the  •'ipole 
initial  transient  motion  and  for  determining  the  effect  of  the  initial  re- 
lease conditions  on  the  steady- state  behavior.  All  of  the  dipole  types 
(test  articles  1 through  13)  were  utilized  for  these  tests. 

6.  Special  Investigations  of  Chaff  Dynamics 

Test  Numbers  400  - 404  were  accomplished  with  a wire-mesh 
disturber  grid  placed  horizontally  across  the  test  chamber.  The  grid 
spacing  was  selected  such  that  there  would  be  a high  probability  of  a 
dipole  striking  the  grid  and  receiving  an  angular  impulse.  The  purpose 
of  this  test  series  was  to  assess  the  dipole  aerodynamic  damping 
moment.  These  tests  were  performed  with  test  article  No.  1. 

Test  Numbers  404  - 415  were  for  the  purpose  of  investigating 
the  effect  of  longitudinal  center  of  gravity  offset  on  the  flight  character- 
istics of  a cylindrical  dipole.  The  2 mil  x 2 inch  dipoles  were  used  for 
these  tests  and  were  weighted  at  one  end  by  dipping  them  in  lacquer  to 
a depth  of  either  5 or  20  diameters. 

Test  Numbers  416  - 421  were  for  the  purpose  of  evaluating  the 
effect  of  geometric  and  mass  asymmetry  on  the  flight  behavior  of  test 
articles  12  and  13.  The  dipole  asymmetries  evaluated  included  twists, 
kinks,  diagonal  end  cuts,  and  various  combinations  of  lead  strip  coating. 

Test  Numbers  422  and  423  were  cluster  drops  of  test  articles 
12  and  13  for  the  purpose  of  photographically  documenting  the  chaff 
initial  dispersal  characteristics. 
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SECTION  V 


DATA  REDUCTION 


1.  Photographic  Data  Processing 

Approximately  600  4 in.  x 5 in.  film  negatives  were  processed 
from  drop  tests  of  individual  chaff  dipoles  in  the  special  aerodynamic 
facility.  These  negatives  were  reviewed  for  image  quality  and  from  them 
175  larger  size  11  in.  x 14  in.  prints  were  produced  for  photogrammetric 
analysis  of  dipcle  position  and  attitude.  In  excess  of  2400  separate  image 
measurements  were  taken  from  these  prints  and  recorded  on  computer 
cards.  Print  reading  was  accomplished  manually  with  the  aid  of  an 
optical  magnifier.  Dipole  position  was  read  with  a metal  scale  to  an 
accuracy  of  approximately  . 005  inches.  Dipole  orientation  was  read  by 
use  of  a drafting  machine  head,  with  an  estimated  reading  accuracy  of 
0.  5 degrees. 

Figures  8 and  9 depict  typical  chaff  dipole  multi-image  photo- 
graphs for  a cylindrical  glass-type  dipole  and  an  aluminum  foil-type 
dipole,  respectively.  Figure  10  shows  multi-image  photographs  of  a 
0.  040-inch  width  foil -type  dipole  taken  from  camera  positions  1 and  3. 
Autorotation  of  the  dipole  about  its  longitudinal  axis  is  clearly  visible  in 
the  latter  photographs. 

2.  Dipole  Trajectory  and  Motion  Data 

A total  of  81  dipole  flight  trajectories  were  converted  into  numer- 
ical time-position-attitude  histories.  This  was  accomplished  using  the 
photogrammetric  procedures  and  data  reduction  computer  program  de- 
scribed in  Appendix  A. 

3.  Aerodynamic  Data  Reduction 

For  each  dipole  trajectory  the  fol'owing  variables  were  deter- 
mined from  numerical  data  processing  of  the  time-position-attitude 
histories: 

• I 

X I 

• I 

Y | velocity  with  respect  to  inertial  reference  system 

• I 

Z | 

Preceding  page  blank 
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p 

q 

r 


dipole  angular  velocity  referenced  to  dipole  moving 
(fixed-plane)  axes 


P 

q 

r 


dipole  angular  acceleration 


CA 

Cn 

CSF 

CM 

CSM 

c* 


dipole  axial  force  coefficient 
dipole  normal  force  coefficient 
dipole  side  force  coefficient 
dipole  pitching  normal  coefficient 
dipole  side  moment  coefficient 
dipole  rolling  moment  coefficient 


a total  angle  of  attack 

4>  angle  of  attack  plane  angular  rotation  parameters 
¥ flight  heading  angle 

y flight  path  descent  angle 

U total  velocity 


Aerodynamic  data  reduction  was  accomplished  using  the  rela- 
tionships described  in  Appendix  B,  and  a Fortran  computer  program 
adapted  to  a CDC  6600  data  processing  machine. 

All  velocities  and  accelerations  were  determined  by  numerical 
moving  polynomial  smoothing  and  differentiation  formulas.  The  number 
of  data  points  for  each  fitted  arc  could  be  selected  by  the  analyst;  either 
a five  or  seven  point  fit  was  employed  in  all  cases.  The  aerodynamic 
coefficients  were  determined  by  direct  solution  of  the  equations  of 
motion  at  the  midpoint  of  each  fitted  arc. 
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SECTION  VI 


RESULTS 


1.  A Qualitative  View  of  Chaff  Dipole  Dynamics 

Both  the  multi-image  photographs  and  the  motion  parameters 
computed  from  the  aerodynamic  data  reduction  program  (see  Section 
IV- 3 and  Appendix  B)  provided  a qualitative  view  of  chaff  dipole  flight 
dynamics. 

The  flight  characteristics  of  a chaff  dipole  were  found  to  be  sig- 
nificantly affected  by  the  cross-sectional  dimensions  of  the  element. 

Slender  Dipoles  Dipoles  with  a diameter  or  width  of  0.008 
inches  or  less  (i.  e. , test  articles  1 through  11)  exhibit  a single  char- 
acteristic motion  in  all  cases.  This  motion  can  be  characterized  as  a 
quasi- steady-state  spiral,  where  the  spiral  angular  rate,  R , and  the 
dipole  azimuthal  angular  rate,  $ , are  nearly  identical  and  either  con- 

stant or  elowly  varying.  The  rate  of  spiral  for  dipole  test  configurations 
1 through  11  varied  from  near  zero  to  about  22  rad/sec.  The  spiral 
rate  is  observed  to  increase  with  either  a decrease  in  dipole  length  or 
an  increase  in  dipole  width  (or  diameter).  Figure  11  shows  frequency 
histograms  of  spiral  rate  for  various  dipole  classificauons. 

The  spiral  is  also  characterized  by  the  rate  of  descent  and  glide 
angle,  Y , as  depicted  in  the  following  sketch. 
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Characterization  of  Chaff  Dipole  Spiral 


The  flight  path  angle  is  closely  related  to  the  trim  angle  of 
attack,  as  defined  in  the  sketch  below. 


dipole 

longitudinal 

axis 


horizontal 


While  the  trim  angle  of  attack  varies  from  neai  zero  to  ninety  degrees, 
the  flight  path  angle  has  a maximum  value  of  the  order  of  20  degrees. 

Both  the  spira’  rate  and  trim  angle  of  attack  are  the  result  of  configura- 
tional asymmetries,  as  described  and  discussed  in  Appendix  C. 

The  final  distinguishing  characteristic  of  the  spiral  motion  (for 
dipoles  with  c S 0.008  inches)  is  that  the  plane  constituted  by  the  dipole 
longitudinal  axis  and  its  velocity  vector  is  always  close  to  vertical.  Thus 
if  the  spiral  were  to  be  compared  to  that  of  a diving  aircraft  in  a highly 
banked  turn,  the  orientation  of  the  dipole  longitudinal  axis  would  corre- 
spond to  the  aircraft  fuselage.  In  other  words,  the  dipole  does  not  ex- 
hibit a "bank  angle"  per  se  but  has  the  appearance  in  flight  of  a gliding 
spear  (see  Figure  8). 

Wide  Dipoles  Dipoles  with  a width,  c , of  0.040  inches  (teat 
articles  12  and  13)  are  dynamically  more  active  than  the  slender  dipoles. 
This  is  apparently  due  to  the  presence  of  a different  flow  regime,  (assoc- 
iated with  the  larger  cross-flow  Reynolds  number)  wherein  vortices  can 
be  shed  from  the  edges  of  the  dipole.  (A  detailed  discussion  of  this 
phenomenon  is  presented  in  Appendix  C).  In  addition  to  the  spiral,  the 
wide  dipoles  can  exhibit  independently  three  types  of  autorotative  motion: 

1)  a large  angle  of  attack  magnus- rotor -type  motion  with  autorotation 
about  the  dipole  longitudinal  axis,  2)  a flat  spin  (or  coning  motion) 
with  autorotation  about  a transverse  axis,  and  3)  a projectile-like 
motion  with  spin  about  the  longitudinal  axis.  The  first  of  these  three 
autorotative  motions  is  shown  in  Figure  10,  while  motions  2)  and  3) 
are  depicted  in  the  photographs  of  Figure  12,  following. 

The  autorotation  rates  of  those  dipoleB  which  experienced  magnus 
rotor  type  motions  were  determined  by  a graphical  motion  fitting  process. 
In  all  cases  the  angular  velocity  was  found  to  be  in  the  range  of  40  - 60 
radians  per  second.  These  rates  correspond  to  values  of  the  nondimen- 
sional  frequency  parameter  u)c/2V  between  0.03  and  0.04. 


2.  Quantitative  Analysis  of  Dipole  Dynamics 

Figures  13-15  present  typical  trajectory  and  motion  history 
data  as  determined  from  the  aerodynamic  data  reduction  program.  Plots 
such  as  these  were  prepared  for  each  of  the  81  dipole  trajectories  which 
were  processed. 

Figure  13  illustrates  the  horizontal  projection  (X-Y  coordinates) 
of  a typical  spiral.  Figures  14  and  15  show  the  time  histories  of  the 
motion  parameters  and  also  the  variation  of  the  aerodynamic  force  and 
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Figure  11  • Multi -Image  Photography  of  . 040-Inch 
Width  Dipoles  in  Autorotative  Motion 
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Figure  14.  Time  Histories  of  Motion  Parameters,  Aerodynamic 
Side  Force  Coefficient,  and  Aerodynamic  Moment 
Coefficients  for  a Typical  Dipole 
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moment  coefficients.  Zero  time  represents  the  initial  image  used  in  the 
data  reduction  and  is  not  the  time  of  release.  In  this  example  the  motion 
parameters  are  slowly  varying,  but  the  motion  does  not  represent  the 
initial  transient  because  the  total  velocity,  U<j>,  is  quasi-steady  at  the 
beginning  of  the  motion  record.  The  correlation  between  the  motion 
parameters  and  also  between  the  motion  parameters  and  the  aerodynamic 
force  coefficients  is  seen  to  be  very  good.  For  example,  as  the  angle- 
of-attack  decreases  there  is  a corresponding  increase  in  the  velocity  and 
a reduction  in  the  normal  force. 

While  the  data  for  most  dipole  descents  displayed  only  a small 
change  in  the  motion  parameters,  the  data  which  have  been  used  for 
illustration  show  that  a single  test  can  provide  usable  conditions  over  a 
range  of  values  of  the  motion  parameters. 

Correlation  of  Descent  Angle  and  Angle  of  Attack  One  of  the 
most  significant  motion  parameter  correlations  is  that  between  the  dipole 
total  angle  of  attack,  a , and  the  glide  angle,  Y . The  results  are  shown 
in  Figure  16,  and  include  data  from  each  of  the  81  dipole  flights  which 
were  analyzed.  The  experimental  results  are  compared  with  the  theoret- 
ical solution  for  the  descent  path  of  a needle-like  body  in  creeping  flow 
(Equation  (C-10)  of  Appendix  C).  It  is  seen  that  the  experimental  data 
closely  agree  wiih  the  theory  both  in  magnitude  and  trend.  The  only 
exception  is  where  rnagnus- rotor -type  motions  were  observed.  Particu- 
larly good  agreement  is  noted  between  the  theory  and  the  experimental 
results  for  the  cylindrical  glass-type  dipoles,  both  of  which  show  a maxi- 
mum glide  angle  of  70  degrees  from  horizontal  at  an  angle  of  attack  of 
35  degrees. 

More  shallow  glide  angles  (i.  e.  , greater  deflections  from  vertical) 
occur  with  the  foil-type  dipoles,  because  these  have  a larger  lifting  force 
due  to  their  greater  projected  surface  area,  but  at  most  the  glide  capa- 
bility of  the  foil-type  dipole  is  only  50  percent  greater  than  that  of  the 
cylindrical  dipole. 

The  above  results  show  that  the  dispersion  characteristics  of 
chaff  dipoles  are  essentially  independent  of  the  dipole  length  and  depend 
only  to  a moderate  extent  upon  the  dipole  cross  section. 

3.  Aerodynamic  Coefficient  Correlations 

The  chaff  aerodynamic  force  and  moment  coefficients  are  defined 
in  Figure  17.  The  coefficients  are  related  to  a right  hand  xyz  body- 
fixed  non- rolling  coordinate  system  which  has  its  origin  at  the  dipole 
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Figure  15.  Correlation  of  Descent  Angle,  y,  with  Angle  of  Attack 


centroid  and  where  the  x axis  coincides  with  the  dipole  longitudinal 
axis,  the  z axis  is  oriented  such  that  xz  plane  always  contains  the 
total  aerodynamic  velocity  vector,  and  the  y axis  completes  the  triad. 
The  normal  force  coefficient,  Cn  , is  perpendicular  to  the  dipole  axis 
and  always  in  the  angle  of  attack  plane,  while  the  side  force  coefficient, 
Csf  > i*  always  perpendicular  to  the  angle  of  attack  plane.  Similarly, 
the  overturning  (or  pitching)  moment  coefficient,  C^ji  corresponds  to 
a rotation  in  the  angle  of  attack  plane,  while  the  side  moment  coefficient, 
Cwm»  relates  to  rotations  perpendicular  to  the  angle  of  attack  plane. 

The  dependence  of  the  aerodynamic  coefficients  on  the  dipole  angular 
orientation  about  its  longitudinal  axis,  which  is  denoted  by  the  angle,  ♦ , 
is  not  considered  directly  in  the  aerodynamic  data  reduction,  because 
this  angle  cannot  be  determined  from  the  photographic  images.  However, 
the  effect  of  is  determined  indirectly  by  the  magnitude  of  Cgp  » since 
this  coefficient  is  zero  for  $ = 0,  and  should  have  a nearly  linear  depen- 
dence on  $ for  < 20  degrees.  The  theoretical  effect  of  4>  on  the  aero- 
dynamic coefficient  is  discussed  in  Appendix  C,  and  is  shown  to  be 
diminishingly  smalt  as  the  dipole  width  decreases. 

Drag  Coefficient  The  relationship  between  the  drag  coefficient 
and  the  previously  defined  normal  and  axial  force  coefficients  is: 

Cj)  - CA  cos  a + Cjsj  sin  a (2) 


Therefore,  for  large  angles  of  attack  (dipole  approximately  perpendicu- 
lar to  the  direction  of  flight) 

Cp  = ' a 90  degrees  (3) 

Values  of  corresponding  to  a > 80  degrees  are  plotted  in  Figure  18 
as  a function  of  a nondimensional  Reynolds  number  parameter 


where 


Rn 

P 

y 

u 


p c u 
y 

air  density 

viscosity 

velocity 

characteristic  dipole  width 
( d if  cylinder;  c if  strip) 


(4) 
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Correlation  of  Normal  Force  Coefficient  with  Reynolds  Number 
for  Large  Angle  of  Attack  Motion  (a  > 80  deg) 


Inspection  of  Figure  18  shows  that  Cjq  correlates  extremely 
well  with  Rjsj  for  the  complete  range  of  test  Reynolds  numbers,  extend- 
ing from  Rn  “ 0.3  to  Rjq  3 50.  The  drag  coefficient  increases  rapidly 
with  decreasing  Reynolds  number,  and  for  the  1-mil  glass  chaff  a maxi- 
mum drag  coefficient  of  36  is  indicated.  The  Reynolds  number  correla- 
tion is  essentially  independent  of  both  the  dipole  length  and  cross-sectional 
shape. 


For  comparison,  theoretical  low-speed  drag  coefficients  for  a 
cylinder  and  plate,  based  on  Oseens  approximations  (Equations  (C-6) 
and  (C-19)  of  Appendix  C)  are  also  shown  in  Figure  18,  along  with  a 
previous  correlation  of  experimental  drag  coefficients  for  cylinders. 
Unfortunately  no  experimental  data  for  flat  plates  could  be  located  in  the 
literature.  The  agreement  between  the  dipole  data  and  these  data  are 
quite  good.  The  only  exception  is  where  the  dipole  is  indicated  to  be  in 
autorotation. 


At  small  angles  of  attack,  where  the  dipole  longitudinal  axis  is 
approximately  aligned  with  the  flow,  the  drag  coefficient  and  axial  force 
coefficients  are  approximately  equal. 

Cd  “ C.A  *»  a ■*  0 (5) 

Values  of  corresponding  to  a 5 30  degrees  are  correlated  with 
Reynolds  number  in  Figure  19.  For  this  correlation  the  characteristic 
length  for  foil-type  dipoles  is  assumed  to  be 


d 


eq 
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TT 


(6) 


Again  dipole  length  and  width  have  practically  no  affect  on  the  correla- 
tion. For  comparison,  the  axial  drag  coefficient  was  computed  from  the 
theory  of  Glauert  and  Lighthill,  as  described  in  Appendix  C.  These 
results,  which  are  also  shown  in  Figure  19,  are  in  good  agreement  with 
the  dipole  drag  data. 

The  present  dipole  aerodynamic  data  apparently  represent  the 
first  experimental  drag  measurements  for  a slender  body  in  axial  flow 
at  very  low  Reynolds  numbers. 

As  a means  of  comparing  the  drag  coefficients  for  axial  and 
normal  type  flow,  the  approximate  fit  to  the  experimental  data  in  Figure 
19  is  re-plotted  in  Figure  18.  It  is  informative  that  as  the  Reynolds 
number  decreases  the  drag  coefficient  is  less  sensitive  to  the  dipole 
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Figure  18.  Correlation  of  Axial  Force  Coefficient  with  Reynolds  Number 
for  Small  Angle  of  Attack  Motions  (a  < 30  deg) 


orientation  and  at  Rj^j  a 0.1  the  drag  coefficients  for  axial  and  normal 
flow  are  nearly  equal.  On  the  other  hand,  at  Rj^j  - 100,  the  drag  coeffi- 
cients for  axial  and  normal  flow  differ  by  more  than  an  order  of  magni- 
tude. 


Normal  Force  Coefficient  at  a Function  of  a The  correlation 
of  the  normal  force  and  social  force  coefficients  with  angle  of  attack  is 
more  difficult,  because  of  the  very  large  influence  of  Reynolds  number. 
For  example,  two  dipoles  of  identical  configuration  at  the  same  angle  of 
attack  can  be  expected  to  have  different  values  of  Cjq  if  their  Reynolds 
numbers  differ.  To  circumvent  this  problem  use  is  made  of  the  results 
from  Stokes  flow  theory  that  the  aerodynamic  forces  are  approximately 
proportional  to  U instead  of  U'-.  Thus,  we  might  expect  for  a given 
value  of  a 


j p S « const  « U (7) 

or  Cft  • U = constant 

This  is  equivalent  to  saying  that  the  normal  force  coefficient  correlates 
as  the  inverse  of  the  Reynolds  number.  If  the  values  of  Cfj  • U for  each 
dipole  cross-section  configuration  are  correlated  separately  with  angle 
of  attack,  the  corresponding  variation  in  cross-flow  Reynolds  number 
for  each  dipole  is  small  and  the  possible  error  due  to  the  Stokes  approxi- 
mation is  also  small. 


Figures  20  and  21  show  the  Cj\j  »U  vs  o correlations  for  the 
cylindrical  glass-type  dipoles  and  foil-type  dipoles,  respectively.  A 
surprisingly  good  correlation  of  Cjsj  • U with  a is  achieved  for  each 
characteristic  dipole  diameter  or  width,  and  the  correlations  are  again 
independent  of  the  dipole  length.  Two  separate  correlations  are  used 
for  the  0.040-inch  width  dipoles,  corresponding  to  whether  autorotation 
about  the  longitudinal  axis  existed  or  did  not  exist.  Through  each  set  of 
Cn*U  values  a curve  proportional  to  sin  a has  been  drawn  and  this 
relationship  is  seen  to  provi<!  a good  fit. 


The  correlation  of  Cjvj  • U with  sin  o constitutes  a validation  of 
cross-flow  theory  for  slender  dipoles,  since  this  is  precisely  the  result 
which  is  obtained  when  the  cross  force  is  equated  to  the  cross  flow 
p d U sin  q 


Reynolds  number 


and  the  local  normal  force  coefficient  is 


assumed  to  be  inversely  proportional  to  the  cross  flow  Reynolds  number. 


Axial  and  Side  Force  Coefficients  Following  the  previous 
approach,  a unique  correlation  of  • U with  a was  also  sought. 
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Results  for  test  articles  9 and  10  are  shown  in  Figure  22.  The  data 
display  a functional  relation  to  cos  a , which  is  again  consistent  with 
cross -flow  theory. 

The  side  force  coefficient,  CgF  * *B  present  only  when  the  dipole 
angular  orientation,  ♦ , is  unsymmetric  with  respect  to  the  cross  flow, 
i.  e.  , HO,  tt/2,  it  ,.. . For  unsymmetric  orientations  the  side  force 
and  normal  force  coefficients  are  related,  theoretically.  Using  equa- 
tions C-20  and  C-21  of  Appendix  C it  can  be  shown  that 


tz  S - 1 

\ 16  TT 

/ _L_\ 

\ 4S2 

/ + Rn 

Us2/ 

mm  “ 


sin  2 4> 


where  the  notation  is  that  of  Appendix  C. 

In  Figure  23  Equation  (8)  is  compared  with  experimental  values 
of  Cgp  and  for  three  representative  foil-type  dipoles  of  identical 
cross-section  (test  articles  9,  10,  and  11).  Inspection  shows  that  a 
majority  of  the  experimental  measurements  correspond  to  values  of  4> 
less  than  20  degrees  and  that  nearly  all  of  the  experimental  values  are 
less  than  the  theoretical  boundary  curve  for  4 = 45  degrees.  The  result 
that  4 is  always  less  than  45  degrees  appears  reasonable  in  view  of  the 
stability  theory  for  ellipsoids  in  an  ideal  fluid,  which  states  that  the  body 
is  stable  only  if  the  motion  is  in  the  direction  of  the  least  axis  (for 
further  discussion  see  Appendix  C). 

The  side  force  coefficient  has  a different  interpretation  for  those 
dipoles  which  experience  autorotation  about  the  longitudinal  axis  and  fly 
at  large  angles  of  attack  (Magnus  rotor  motion).  In  such  cases,  the  side 
force  coefficient  corresponds  to  the  classical  magnus  force,  in  accord- 
ance with  the  following  sketch. 


WI  1^1W¥1W»1  V.  — . 

with  Angle  of  Attack  (Foil-Type  Dipoles) 


X 
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For  all  of  those  dipole  flights  which  experienced  magnus  rotor 
motions  the  side  force  coefficient  had  a magnitude  of  about  0.25.  Defin- 
ing the  magnus  force  coefficient  as 


3 C 


SF 


0)  c 
2 U 


(9) 


values  of  Cjsjp  from  about  6.  0 to  9.  0 were  computed.  No  explanation 
can  be  given  for  these  large  values  of  , which  exceed  in  magnitude 

experimental  measurements  for  rotating  cylinders. 


Aerodynamic  Moment  Coefficients  At  very  low  Reynolds  num- 
bers theory  indicates  that  the  moment  due  to  fluid  pressure  vanish  if  the 
body  has  three  mutually  perpendicular  planes  of  symmetry.  Direct 
measurements  of  the  dipole  moment  coefficients  C\i  and  Cgj^j  confirm 
the  extremely  small  rr^pmtude  of  the  aerodynamic  moments  (see 
Figures  14  and  15).  Howev.”  two  significant  sources  of  aerodynamic 
moment  .emain;  1)  the  moments  du“  to  configurational  asymmetries, 
i.  e.  , bending,  twist,  etc.,  a:  * 2)  the  viscous  moments  caused  by 
dipole  angular  motion.  Since  the  measured  moment  coefficients  reflect 
all  three  moment  contributions,  a separate  determination  of  each 
moment  contribution  is,  in  general,  impossible.  Consequently,  indirect 
methods  of  assessing  the  probable  magnitude  of  the  aerodynamic  moment 
contributions  were  employed.  The  ba^ic  approach  involved  calculation 
of  the  viscous  damping  moment,  and  then  inference  of  the  static  moment 
contributions  of  the  symmetric  and  asymmetric  dipole.  The  method  of 
estimating  the  viscous  damping  is  described  in  Appendix  C for  dipole 
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rotations  about  either  the  y or  z axes  of  Figure  17.  Using  estimated 
damping  moments  the  dipole  equations  of  motion  were  then  re-evaluated 
and  the  remaining  static  moment  coefficients  determined.  Figure  24 
illustrates  the  values  of  Cm  and  Cgj^  determined  in  this  fashion,  using 
motion  data  for  a large  selection  of  dipole  configurations. 

The  following  tentative  conclusions  are  drawn  from  this  data: 

1)  The  static  moment  coefficient,  Cm  , for  a symmetric  dipole 
is  extremely  small  and  probably  not  larger  than  0.005  for 
glass-type  dipoles  or  larger  than  0.015  for  foil-type  dipoles. 

2)  The  side  moment  coefficient  CgM  is  significantly  larger 
than  the  overturning  moment  coefficient  Cm  . 

3)  For  foil-type  dipoles  both  Cm  and  CgM  attain  maximum 
values  at  an  intermediate  angle  of  attack  in  the  neighborhood 
of  a = 45°,  but  for  glass-type  dipoles  CgM  is  maximum 

at  small  angles  of  attack. 

The  larger  magnitude  of  CgM  compared  to  Cm  for  the  foil  type  dipoles, 
is  probably  due  to  the  presence  of  twist,  which  has  a large  contribution 
to  CgM  but  only  a small  contribution  to  Cm  • 

As  a further  means  of  investigating  the  magnitude  of  the  static 
overturning  moment,  several  drop  tests  of  test  article  1 were  accom- 
plished with  intentional  center -of-mass  offset.  The  smallest  offset 
tested  ( Ax  = 0.  01  £)  was  sufficient  to  trim  the  dipole  to  an  angle  of 
attack  approaching  zero.  From  this  it  was  inferred  that  the  maximum 
static  moment  of  the  symmetric  dipole  could  not  have  exceeded  0.03, 
thus  confirming  the  test  data  shown  in  Figure  24. 

( Finally,  a special  experiment  was  conducted  to  provide  a direct 
evaluation  of  the  aerodynamic  damping  moment  in  pitch  (see  Section  IV-6) 
for  comparison  with  Equation  (C-34)  of  Appendix  C.  By  direct  measure- 
ment it  was  found  that  M = -4.7  x 10"^  ft  - lb  and  from  Equation  (C-34)  a 
value  of  M = -2.4  x 10“^  ft-lb  was  obtained.  Thus,  the  damping  moment 
predictions  are  believed  to  be  accurate  within  a factor  of  two,  which  is 
reasonable  considering  the  extremely  small  magnitude  of  the  moments. 

i 

A more  extensive  discussion  of  the  moments  due  to  configura- 
tional asymmetry  can  be  found  in  Appendix  C, 
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Figure  23.  Correlation  of  Static  Overturning  Moment  Coefficient 

and  Static  Side  Moment  Coefficient  with  Angle  of  Attack 
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4. 


Chaff  Cluster  Dispersion  Testa 


An  experimental  investigation  of  chaff  cluster  dispersion  was 
accomplished  as  an  adjunct  to  the  basic  aerodynamic  test  prog:  am.  The 
purpose  of  the  cluster  drops  was,  1)  to  assess  the  affect  of  dipole 
spiraling  over  a time  span  encompassing  many  cycles,  2)  to  determine 
the  relative  dispersion  contributions  of  the  transient  and  steady-state 
descent  phases,  and  3)  to  examine  the  possible  contribution  of  atmo- 
spheric turbulence  to  chaff  dispersion  in  an  airship  hangar  environment. 
The  tests  were  accomplished  using  chaff  cluster  of  either  IOC  or  200 
dipoles.  All  drops  from  2 5 ft.  height  and  above  were  accomplished  in 
the  MCAS  airship  hangar  at  Santa  Ana,  Calif.  To  minimize  air  circu- 
lation the  tests  were  made  between  3-6  A.  M.  with  all  the  hangar  doors 
closed. 


Dispersion  Patterns  Ground  level  dispersion  patterns  for  a 
representative  small  chord  foil-type  dipole  (test  article  No.  9),  as 
recorded  from  cluster  releases  at  25,  50,  and  100  ft.  above  ground 
level  are  shown  in  Figure  25.  Similar  data  were  obtained  for  test 
article  No.  12.  The  recovery  factor  for  these  drops  varied  from  100 
percent  at  the  lowest  drop  height  to  about  90  - 95  percent  for  the  100 
ft.  drop  height.  An  attempt  to  measure  the  dispersion  of  the  glass-type 
dipoles  was  unsuccessful  because  the  recovery  was  less  than  50  percent 
for  the  minimum  drop  height. 

Statistical  analysis  of  the  impact  patterns  was  accomplished  using 
the  cumulative  frequency  distribution  for  the  radial  deflection  from  the 
mean  center  of  impact.  The  radial  deflection  statistics  were  found  to 
closely  match  the  two-dimensional  normal  distribution  for  cumulative 
frequencies  up  to  about  80  percent,  but  for  large  dispersions  the  data 
showed  non-Oussian  trends.  Values  of  the  probable  dispersion  (50th 
percentile)  are  shown  as  a function  of  release  height  in  Figure  26  for 
both  wide  and  narrow  types  of  foil  chaff. 

The  test  results  indicate  that  for  release  heights  above  approxi- 
mately 25  ft.  the  dipole  dispersion  is  essentially  constant,  if  the  drops 
are  accomplished  in  a protected  environment.  Below  25  ft.  , the  disper- 
sion increases  rapidly  with  increasing  release  height.  For  a drop  height 
of  100  ft.  there  is  no  significant  difference  in  dispersion  between  the 
two  dipole  configurations  tested. 

Spiral  Analysis  Using  probable  values  of  spiral  rate,  descent 
velocity  and  flight  path  angle  (as  determined  from  the  instrumented  drop 
tests),  the  probable  values  of  turning  radius,  R,  and  probable  dispersion. 
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b)  Release  Height  B 50  Feet 
Figure  24.  (Continued) 
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c)  Release  Height  a 100  Feet 
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d)  Release  Height  * 100  Feet  (Repeat) 


Figure  24.  (Continued) 
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Large-Chord  Foil  Dipoles  Small-Chord  Foil  Dipoles 
O 100  element  cluster  □ i00  element  cluster 
Cf  200  element  cluster  rf  200  element  cluster 


Probable  Radial  Dispersion  from  Mean  Center  of  Impact  as  a Function 
of  Cluster  Release  Height  (Foil-Type  Dipoles) 


r , for  a steady-state  spiral  were  computed  from  Equations  (C-42)  ano 
(C-47)  of  Appendix  C.  The  results  are  given  below: 


Test 

Dipole  Dimensions 

R 

r 

Article 

inches 

ft. 

ft. 

9 

. 006  x . 00045  x 1. 78 

0.  21 

0.  27 

12 

. 040  x . 00045  x 1. 12 

0.  32 

0.41 

The  probable  dispersion  due  to  steady- state  spiralling  as  given 
by  the  table  above  is  seen  to  be  significantly  less  than  the  observed  dis- 
persion data  shown  in  Figure  25.  It  can  therefore  be  concluded  that 
i rider  quiescent  atmosphere  conditions  the  release  transient  is  the 
primary  source  of  dipole  dispersion  for  descent  distances  of  the  order 
of  100  feet. 
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SECTION  VII 


SIX-DEGREES-OF-FREEDOM  MOTION  SIMULATION 

1.  6 -DOF  Computer  Program 

A specially  modified  6-DCF  computer  program  was  prepared  for 
simulation  of  the  complete  motion  of  single  or  multiple  dipoles. 

The  basic  computer  program  from  which  the  final  program  was 
derived  is  the  Alpha  Research  "Extended  Capability  Magnus  Rotor  and 
Ballistic  Body  6-DCF  Trajectory  Program"  which  is  documented  in 
References  2,  3,  and  4.  The  basic  program  has  such  features  as 

• all  attitude  motion  prediction 

• self-adjusting  integration  schemes 

• option  for  either  body- fixed  or  fixed-plane  axes 

• three  parameter  aerodynamic  coefficient  tables  plus 
aerodynamic  dependency  upon  roll  orientation 

• provision  for  aerodynamic,  geometric,  and  inertial 
asymmetries 

• Monte  Carlo  operation,  i.  e.  , random  selection  of  initial 
motion  parameters  and  asymmetries 

• provision  for  modeling  of  an  initial  cluster  break-up 

• versatile  input  format  in  NAMELIST  notation. 


^ Alpha  Research,  Inc.  , "User's  Manual:  Extended  Capability 
Magnus  Rotor  and  Ballistic  Body  6-DOF  Trajectory  Program," 
Report  No.  AFATL-TR-70-40,  May  1970. 

^ Alpha  Research,  Inc,,  "Amended  User's  Manual:  Extended 
Capability  Magnus  Rotor  and  Ballistic  Body  6-DOF  Trajectory 
Program,"  Alpha  Research  ReportNo.  71-2,  26  March  1971. 

* Alpha  Research,  Inc.  , "6-DOF  Monte  Carlo  Trajectory  Program", 
Alpha  Research  Report  No.  72-0089-10,  29  September  1972. 
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For  chaff  simulation  the  following  additional  features  were  provided; 


/ 

f 

I 


• atmospheric  wind  and  turbulence  model 

• aerodynamic  coefficient  dependence  upon  Reynolds  number 

• special  output  format 

• improved  integration  controls. 

The  atmospheric  wind  model  has  options  for  either  steady  winds 
of  arbitrary  direction,  or  an  arbitrary  vertical  wind  shear  (of  fixed 
directional  heading)  with  superimposed  isotropic  stochastic  turbulence. 
The  turbulence  model  incorporates  both  Eulerian  and  Lagrangian  scales, 
i.  e. , both  spacial  and  time  dependency.  The  correlation  functions  are 
tailored  to  chaff  descent  conditions. 

A complete  description  of  the  turbulence  modeling  is  presented  in 
Appendix  D.  Selection  of  the  turbulence  option  in  the  6-DOF  trajectory 
program  requires  only  four  additional  parameters,  which  are  input  as 
table  functions  of  altitude.  Recommended  values  of  the  turbulence 
parameters,  as  well  as  a discussion  of  atmospheric  turbulence  measure- 
ments, may  be  found  in  Appendix  E. 


Aerodynamic  Coefficient  Data 


Complete  6-DOF  aerodynamic  data  packages  have  been  prepared 
for  two  representative  dipole  configurations:  the  1 mil  x 1-inch  glass- 
type  dipole  and  . 006  x . 00045  x 1.78-inch  foil-type  dipole. 


The  aerodynamic  coefficients  for  the  1 mil  x 1-inch  dipole  are 
presented  in  Table  3 as  a function  of  both  angle  of  attack  and  Reynolds 
number,  Rjyj  ( £)*.  The  aerodynamic  coefficient  notation  is  consistent 
with  Figure  17.  The  force  coefficients,  Cx  and  Cj$,  as  well  as  the 
damping  derivative  ®»  Cmq  and  Cnr  , are  baeed  on  data  of  Figures  18 

and  19,  and  the  following  relationships; 


* The  6-DOF  trajectory  program  accommodates  only  one  characteristic 
length  and  the  dipole  length,  l , has  been  selected  fc-  this  parameter, 
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Table  3.  BASIC  AEROBALLISTIC  COEFFICIENTS  FOR  CYLINDRICAL  GLASS  DIPOLE  (O.OOl  x 1.0  INCHES) 
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(C*U 


cos  a 


(CN>a= 


n/2 


^CN^7t/2 


,^A/2 


sin  a 
sin^  a 


(ID 

(12) 

(13) 


The  latter  two  formulas  correspond  to  Equations  (C-35)  and  (C-38)  of 
Appendix  C.  The  maximum  value  of  the  overturning  moment  coefficient 
was  taken  to  be  0.005  based  on  the  data  of  Figure  24.  This  value  was 
assumed  to  be  appropriate  for  the  nominal  descent  Reynolds  number. 

The  moment  coefficient  was  further  assumed  to  decrease  with  decreasing 
Reynolds  number  in  direct  proportion  to  the  change  in  Reynolds  number. 
The  roll  damping  coefficient,  Cj^,  was  estimated  using  the  following 

relationship,  which  is  similar  to  that  for  rotating  spheres. 


% ■ 

k ir  d 

Ru  (i,)  a 

where 

d * 

dipole  diameter 

l = 

dipole  length 

Rjj  ( A ) = 

length  Reynolds  number 

(14) 


The  body-fixed  moment  coefficient  is  used  to  represent  the  effect 

of  random  longitudinal  bend.  The  moment  coefficient  values  are  normal- 
ized such  that  the  standard  deviation  of  bend  corresponds  to  a 6 value 
of  unity.  The  corresponding  value  of  the  dipole  longitudinal  bend  param- 
eter, R/a  (as  described  in  Appendix  C,  Section  3)  was  arbitrarily 
established  as  5.0.  The  moment  coefficient  due  to  longitudinal  bend 
includes  both  the  axial  and  normal  force  contributions  as  given  by 
Equations  (C-28)and  (C-31),  respectively. 

The  aerodynamic  coefficients  for  the  0.006  x .00045  x 1.78-inch 
dipole  are  presented  in  Table  4.  In  addition  to  the  functional  dependence 
upon  a and  Rjsj , the  functional  dependence  upon  the  roll  angle,  <S>  , is 
also  incorporated  through  a set  of  additional  harmonic  coefficients 
defined  as  follows: 
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Table  4.  BASIC  AEROBALLISTIC  COEFFICIENTS  FOR  FOIL  DIPOLE  (0.006  x .00045  x 1.78  INCHES) 


1 


NoUs  CNPA  r CM  PA 


Table  4.  ADDITIONAL  BODY-FIXED  AERODYNAMIC  COEFFICIENTS  FOR  FOIL  DIPOLE 
(0.006  x .00045  x 1.78  INCHES)  (Continued) 


Not*!  PHSFD  : PHMD  = PHSM2D  = 90.0;  PHND  - PHSM1D  = PHUD 


ACjj  = [CHl  (a,  Rn)]  sin  2 $ 

(15) 

ACgp  = [Cgp^  (a,  Rjj)]  sin  (2  $ + tt/2 ) 

(16) 

acM  * [CMX  (a»  %)]  sin  (2  $ + tt/2) 

(17) 

ACsm  = [Csm1  (a.  %))  sin  2 * 

(18) 

♦ [Csm2  (a.  Rjj)]  COS  2 4> 

AC£  = lC£#  (a,  Rjj)]  sin  2 $ 

(19) 

The  functional  form  of  the  above  equations  is  derived  from  Appendix  C, 
Section  2,  and  provides  for  not  only  the  rotational  symmetry  of  the  foil- 
type  dipole,  but  also  for  twist  and  longitudinal  bend  effects  as  given 
by  Table  C-l.  The  coefficient  contributions  from  twist  are  normalized 
to  a standard  deviation  e = 10  degrees,  while  the  contributions  from 
longitudinal  bending  are  normalized  to  a bend  ratio  R/a  = 5.0. 

The  determination  of  ihe  coefficient  values,  with  the  effect  of  ♦ 
included,  requires  special  attention.  Consider,  for  example,  Cj*j  and 
Cnj  . In  accordance  with  Table  C-l 

Cjj  (a,  Rjj ) + (a,  Rjj,  4>)  * sin  a + k2  sin  a cos  2 <t 


Therefore 

Cjj  (a,  Rjj)  = kj  sin  a 

(20 

Cjj^  (a,  Rjj,  <J>)  = k2  sin  a cos  2 <J> 

(21) 

or 

Cjj1  (a,  Rjj)  = k2  sin  a 

(22) 

However,  the  experimental  data  represent  the  maximum  normal  force 
( 4>  -*•  0)  and  include  both  the  Cjvj  and  Cj\j.  contributions.  Thus,  for 
a = tt/2 

^N^meesured  = = ^1  + ^2  ^3) 

maximum 

To  obtain  kj  and  k2,  separately,  an  additional  relationship  is  required, 
such  as  the  ratio  of  the  minimum  and  maximum  normal  force.  At  suffi- 
ciently low  Reynolds  number  the  ratio  of  (Cjvj)mjn/(Cjsj)max  can  be 
obtained  with  the  aid  of  Equation  (C-20).  For  small  Reynolds  numbers 
this  ratio  is  found  to  have  a value  of  about  0,8.  Defining  the  ratio 
(cN)min/(cN)max  as  c « is  easily  shown  that 
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(24) 


kl  = <cN>tt/2 

k2  = (CN)ff/2  L~  (25) 

For  convenience,  a summary  of  all  the  k's  and  K's  used  in 
computing  the  aerodynamic  coefficient  tables  is  given  in  Table  5. 
Insufficient  data  were  availat’ * to  determine  K2;  k,|  was  neglected 
because  of  its  small  contribution. 

3.  Coefficient  Tables  for  Alternate  Chaff  Lengths 

With  appropriate  adjustment  of  reference  area,  reference  length, 
and  the  Reynolds  number  argument  values,  coefficient  values  in  Tables 
3 and  4 may  be  used  without  further  modification  for  alternate  chaff 
lengths  provided  that 

a)  the  dipole  slenderness  ratio  remains  large 

b)  the  longitudinal  bend  ratio,  R/a,  remains  constant 

c)  the  twist  angle,  e , remains  constant. 

Defining  a new  chaff  lenrth,  , the  new  Reynolds  number 
argument  values  are 

Rn  (V)  * Rn  (*)  (26) 

The  aerodynamic  reference  area  is  either 

S « dl'  or  S « c*' 

(for  cylindrical  and  strip-type  dipoles,  respectively)  and  the  aerodynamic 
reference  length  is  V . 
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4. 


Monte  Carlo  Simulations 


To  illustrate  the  use  of  6- DOF  simulation  for  chaff  dipole  motion 
analysis,  a Monte  Carlo  trajectory  analysis  was  accomplished  for  the 
. 006  x . 00045  x 1. 78-inch  foil-type  dipole  using  the  modified  computer 
program  described  in  Section  (VI- 1).  For  this  series  of  trajectory 
simulations,  a quiescent  atmosphere  without  wind  or  turbulence  was 
assumed. 

Physical  data  for  the  0.  006  x . 00045  x 1.  78-inch  foil-type  dipole 
were  taken  from  Table  2,  and  aerodynamic  coefficient  data  from  Table  4. 
The  dipoles  were  assumed  to  be  released  from  a fixed  position  with  a 
random  initial  orientation.  The  dipoles  were  further  assumed  to  have 
Gaussian  center-of-mass  statistics,  with  the  standard  deviations  of 
center-of-mass  offset  corresponding  to  values  of  length,  width,  and 
thickness  of  0.0G5JI,  0.01c,  and  0.  02t,  respectively.  The  standard 
deviation  of  dipole  longitudinal  bend  was  taken  as  R/a  = 5.  0 (i.  e.  , 

R = 4.  45  inches)  and  the  standard  deviation  of  element  twist  was 
assumed  to  be  10  degrees. 

The  6-DOF  motion  calculations  were  performed  with  the  computer 
program  set  for  the  Monte  Carlo  mode.  For  this  series  of  simulations 
the  integration  time  interval  was  automatically  adjusted  and  was  nomin- 
ally in  the  range  of  0.001  - 0.004  seconds. 

Figure  27  depicts  the  XY  coordinates  of  the  ten  sample  trajec- 
tories which  were  computed  for  standard  sea  level  conditions.  These 
trajectories  represent  a flight  time  of  4.  0 seconds.  The  projected 
orientation  of  each  dipole  at  a representative  point  along  its  trajectory 
is  illustrated  as  a means  ot  further  portraying  the  computed  motion. 

The  computed  descent  velocity  varied  from  1.  1 to  1.9  ft/sec,  which 
compare  with  measured  velocities  of  1.  1 to  3.4  ft/sec  for  this  dipole 
configuration. 

Although  no  statistical  analyses  of  the  computed  motions  were 
accomplished,  qualitatively  the  data  appear  quite  similar  to  that  ob- 
tained in  the  aerodynamic  test  chamber. 

5.  Simulation  of  Dipole  Flight  in  a Turbulent  Atmosphere 

Several  six-degrees-of-freedom  simulations  of  dipole  flight  under 
turbulent  atmospheric  condition  were  accomplished  with  the  specially 
modified  6-DOF  computer  program,  using  representative  low  and  high 


altitude  turbulence  intensities.  Motion  simulations  were  accomplished 
for  both  the  1 mil  x 1-inch  and  .006  x .00045  x 1.75-inch  dipole  con- 
figurations. As  with  the  Monte  Carlo  simulations,  random  initial 
conditions  and  configurational  asymmetries  were  introduced. 

Figure  28  illustrates  typical  motion  data  for  the  . 006  x . 00045  x 
1.78-inch  foil-type  dipole  in  low  altitude  turbulence  starting  at  an  initial 
height  50  ft.  above  ground  level.  The  turbulence  parameters  from  Table 
E-2  were  used  as  input  to  the  computer  program  for  this  simulation. 

It  is  seen  from  Figure  28  that  turbulence  affects,  primarily,  the 
inertial  velocity  components,  such  as  the  vertical  descent  rate  which  is 
shown.  The  aerodynamic  velocity,  i.  e. , the  velocity  of  the  element 
minus  the  wind  velocity,  fluctuates  only  slightly.  A very  significant 
result  is  that  the  dipole  attitude  variation  (as  represented  by  the  Euler 
angles  0 and  ^ ) is  quite  steady,  even  though  the  orientation  of  the 
flow  relative  to  the  dipole  varies  widely.  Thus,  the  angular  motion  of 
the  dipole  is  practically  unaffected  by  turbulence.  This  is  probably  due 
to  the  aerodynamic  damping,  which  is  extremely  large  in  comparison  to 
forcing  moments.  The  translational  motion  of  the  dipole  with  respect  to 
an  observer  moving  with  the  mean  wind,  on  the  other  hand,  would  appear 
erratic,  and  might  be  termed  a "jitter"  motion. 
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SECTION  VIII 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


! 

I 


The  aerodynamic  characteristics  of  thirteen  distinct  dipole  con- 
figurations, varying  in  construction,  cross-section,  and  length,  were 
determined  from  instrumented  drop  tests  of  single  elements  in  an 
enclosed  test  chamber.  The  dipole  motion  was  recorded  photographic- 
ally with  three  orthogonal  still  cameras  equipped  with  p — >chronized 
rotating  shutters. 

The  dynamic  behavior  of  a chaff  *!  _t  A-as  found  to  depend 
greatly  upon  its  principal  cross-sectional  dimension.  Dipoles  with 
cross-sectional  dimensions  equal  or  less  than  about  0.008  inches 
exhibited  a singular  characteristic  spiral  motion,  even  though  there 
were  large  differences  in  spiral  rate  and  trim  orientation.  In  all  cantts 
the  spiral  rate  and  dipole  yawing  rate  were  essentially  equal.  Dipoles 
with  a characteristic  cross-sectional  dimension  of  0.040  inches  were 
dynamically  more  active  and  sometimes  displayed  auto  rotations  about 
either  the  longitudinal  axis  or  a transverse  axis. 

Aerodynamic  force  coefficients  for  all  of  the  dipole  configura- 
tions could  be  correlated  with  the  cross-Reynolds  number.  With  appro- 
priate choice  of  parameters  both  the  normal  force  and  axial  force 
coefficients  were  found  to  be  independent  of  dipole  length  and  cross- 
sectional  shape,  and  dependent  only  upon  angle  of  attack  and  Reynolds 
number.  Agreement  of  the  dipole  drag  data  with  theory  and  other 
experimental  data  on  cylinders  at  low  Reynolds  number  is  good.  For 
the  smallest  dipole  diameter,  a maximum  drag  coefficient  of  38  was 
determined. 

The  static  aerodynamic  moments  acting  upon  the  dipoles  with 
characteristic  cross-section  dimensions  equal  or  less  than  0.008  inches 
were  found  to  be  extremely  small  and  due  primarily  to  configurational 
asymmetries  such  as  longitudinal  bending  and  twist.  Procedures  for" 
estimating  the  forces  and  moments  due  to  asymmetry,  and  the  angular 
velocity  damping  were  derived,  and  these  methods  were  generally  sub- 
stantiated by  the  test  data. 

Using  both  the  experimental  data  and  aerodynamic  theory,  6- DOF 
aerodynamic  data  packages  were  prepared  for  two  representative  dipole 
configurations.  These  aerodynamic  data  packages  were  used  in  conjunc- 
tion with  a specially  modified  6- DOF  computer  program  to  simulate  the 
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dipole  motion  under  both  quiescent  and  stochastic  turbulent  atmospheric 
conditions.  Both  single  element  and  Monte  Carlo  type  simulations,  with 
random  initial  conditions  and  dipole  configurational  asymmetries,  were 
accomplished.  The  computed  dipole  trajectories  and  motion  are  in 
qualitative  agreement  with  the  experimental  motion  data.  The  effect  of 
turbulence  on  the  dipole  motion  was  found  to  be  primarily  a translational 
jitter.  The  dipole  angle  motion  is  relatively  insensitive  to  turbulent 
wind  fluctuations. 

The  dispersive  char  .eristics  of  chaff  dipole  clusters,  for  drop 
-.u'.  8 - 100  above  ground  level,  were  determined  from  drop 
tests  in  an  airs*..,.-  hangar  under  quiescent  atmospheric  conditions.  For 
this  range  of  d.-op  heights  chaff  dispersion  is  due  largely  to  the  dipole 
initial  motion  transient  and  the  dispersion  is  nearly  constant  for  drop 
heights  between  25  and  100  feet. 

The  present  effort  has  provided  a general  understanding  of  chaff 
dipole  aerodynamics  and  dynamic  behavior.  It  has  also  provided  the 
mjans  for  analytic  simulation  of  chaff  dipole  motion  in  the  presence  of 
stochastic  environments.  However,  much  remains  to  be  done  in  terms 
of  assessing  the  relationship  between  chaff  element  dynamics  and  chaff 
cloud  electromagnetic  performance.  Specific  areas  of  work  or  studies 
which  should  be  accomplished  include: 

1)  Monte  Carlo  simulations  with  large  sample  size 

2)  more  detailed  study  of  the  relationship  between  dipole 
asymmetry  and  its  motion 

3)  investigation  of  simplified  aerodynamic  models  as  a 
means  of  reducing  computer  time  requirements  for 
6-DOF  simulations 

4)  evolvement  of  improved  dipole  configurations  on  the 
basis  of  analytical  simulation  of  chaff  cluster  dynamic 
behavior  and  electromagnetic  performance. 
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APPENDIX  A 


1 

I 

PHOTOGRAMMETRIC  DATA  REDUCTION  j 

* 

The  descent  trajectories  of  the  aerodynamic  test  dipoles  are  • 

recorded  by  three  still  cameras  with  rotating  shutters  (Figure  A-l).  j 

These  cameras  are  so  arranged  as  to  provide  orthogonal  views  of  the  \ 

drop  test  chamber  from  each  of  two  camera  pairs  (cameras  1 and  2 or  j 

cameras  1 and  3).  The  instantaneous  position  and  orientation  of  the  I 

descending  dipole  is  documented  by  multi-image  tracings  on  each  film  j 

negative.  Through  photogrammetric  analysis  the  true  pitch  and  yaw  < 

attitudes  of  the  dipole  are  described  in  standard  Euler  spherical  angles  1 

and  the  dipole  space  position  is  determined  in  rectangular  space  coordi-  ' 

nates.  The  photogrammetric  data  reduction  procedure  also  removes  all  j 

optical  parallax  errors.  ! 

{ 

Because  the  dipole  position  and  attitude  can  be  determined  by  - 

either  of  two  camera  pairs,  the  above  described  motion  recording  system 
has  a valuable  redundancy.  This  redundancy  feature  can  be  exploited  1 

when  the  image  quality  is  poor  from  one  camera  station,  or  when  the  \ 

aspect  of  one  camera  is  unfavorable  for  accurate  attitude  and  position  ! 

determination. 

1.  Determination  of  Dipole  Position  and  Attitude  (Cameras  1 and  3)  | 

Cameras  1 and  3 comprise  one  of  the  orthogonal  sets,  with  the  j 

camera-3  axis  positioned  such  that  it  is  0O  degrees  off  vertical.  The 

camera-1  axis  is  aligned  horizontally.  The  determination  of  dipole  posi-  Ij 

tion  and  attitude  from  cameras  1 and  3 is  slightly  more  involved  than  I; 

with  cameras  l and  2.  ; ; 

Let  the  location  of  the  lens  of  camera  3 be  defined  as  X = b, 

Y = d,  and  Z = e.  The  XY  plane  bisects  the  center  of  the  4-ft  cube  test 
chamber,  and  the  lens  axes  of  cameras  1 and  2 both  lie  in  this  plane.  > 

The  lens  of  camera  1 is  located  at  X = a,  Y = b,  and  Z = 0 in  accord- 
ance with  Figure  A-l. 

Let  the  chaff  dipole  be  located  in  the  4-ft-cube  drop-test  zone  at 
Xp,  Ypi  Zp,  and  let  the  cameras  be  located  as  indicated  previously.  A 
line  through  Camera  1 and  the  geometric  center  of  the  particle  has  the 
equation: 
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Figure  A-l.  Coordinate  System  for  Motion  Recording  Cameras 
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X-a 


Y-b 


Z-0 


(A- 1 ) 


\ 

\ 


X 


Xp-a  " Yp-b  " Zp-0 


The  intersection  of  this  line  and  the  YZ  pla..*e  «a  at: 
Y H \ , Z = Clf  X = 0 


Thus 


(A-2) 


(A-3) 


Consider  a line  from  camera  3 to  the  particle  as  having  the  equation: 


X-b  Y-d  Z-e 
Xp-b  = Yp-d  = Zp-e 


(A-4) 


A new  coordinate  system  is  now  defined  *n  which  plane  X3Y3  is  perpen- 
dicular to  the  center-view  axis  of  camera  3. 


V 
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The  transformation  matrix  is: 


X 

’10  0 

V 

Y 

= 

0 c22  c23 

Y3 

Z 

.0  c32  c33_ 

-Z3. 

Substituting  these  values  into  Equation  (A-4)  results  in 

X3-b  C22Y3  + ^23^3  " ^ C32Y3  + C33Z3  " e 


Xp-b 


Yp-d 


zp-e 


(A- 5) 


The  image  in  the  X3Y3  plane  when  looking  through  camera  3 
appears  to  be  at 

X3  = 53.  Y3  = n3  , Z3  = 0 

Substituting  these  values  to  Equation  (A-5)  and  combining  with  Equation 
(A-2)  Xp,  Yp,  and  Zp  are  determined  as  follows: 


rij  - d + b/f3 


Yp  = b - (Xp-a)  ( nrb)/a 


(A-6) 


(A-7) 


Cl 

T (Xp-a) 


(A-8) 


Equations  (A-6),  (A-7),  and  (A-8)  are  used  to  locate  the  position  of  the 
dipole. 


Let  0 and  <p  be  the  Euler  angles  describing  the  attitude  of  the 
long  axis  of  the  dipole.  Also  define  03  as  the  angle  seen  effectively 
between  the  filament  axis  and  the  Y axis  in  the  YZ  plane  in  accordance 
with  the  sketches  on  the  following  page. 
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Angular  Coordinates 


The  following  direction  vectors  are  now  defined: 

Ip  = direction  vector  through  the  filament  axis 

-k.  -+  -► 

ip  r sin  0 co 8 $ i + sin  0 sin  <J>  j + cos  0 k 
itj  = direction  vector  from  the  particle  to  camera  1 
= (Xp-a)  i + (Yp-b)  J + Zpk 

If  the  coefficients  of  a normal  vector  are  n = x , then  the  plane 
which  contains  li  and  f.p  can  be  determined.  The  intersection  between 
this  plane  and  the  YZ  plane  {X  = 0)  is  a line  with  the  slope: 


m = tan  = 


(Xp-a)  cos  0 - Zp  sin  0 cos 
(Xp-a)  sin  0 sin  $ - (Yp-b)  sin  0 cos  $ 


(A-9) 


Rearranging, 

“p-** 

tan  [(Xp-a)  tan  sin  $ -(Yp-b)  tan  a]_  cos  <J>  + Zp  cos  $ ] 


Xp-a 


Now  consider  the  X3Y3Z3  grid  space  where 
angle  between  the  dipole  axis  and  the  X3  axis  in  the 


013  is  defined  as  the 
X3Y3  plane. 
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In  the  X3Y3Z3  coordinates,  Camera  3 is  located  at  (b3g3h3),  where: 

b3  = b 

83  = c22d  + c32e 
h3  = C23d  + C33e 

£3  r vector  from  the  particle  to  Camera  3 

£3  = (Xp3  - b3)  i3  + (Yp3  - g3)  J3  + ( Zp3  - h3)  k3 


The  direction  vector  £p  was  given  previously,  but  transformed  to 
X3Y3Z3  coordinates  becomes: 


£p  = sin  6p  cos  <Jp  i3  + (C22  sin  OpSin^  +C32cos  0p)  J3 


+ (C23  sin  0p  sin  $r  + C33  cos  0p  ) k3 
Thus,  the  normal  is  given  by: 


J3 


XP3  "b3 


YP,-83 


ZP3'h3 


sin  0 cos  p 


C22  sin  6 sin  Q 

+ C32  cos  0 


C23  sin  0 sin  <f> 


+ C33  cos  0 


By  similarly  using  the  coefficients  of  n , a plane  can  be  defined 
perpendicular  to  n . This  plane  intersects  the  X3Y3  plane  at  Z3  = 0 . 
The  slope  of  that  line  is: 


mi  = tan  a,  = - 


coeff  ( i3  term) 


coeff  ( term) 


■M 


1 
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tan  a- 


Yp4  (d2C23  + d3C33)  - Zp4  (d2C22  + d3C32) 
Xp4(d2C23  + d3C33)  - Zp4dt 


where 


dj  = sin  e cos  $ 
d2  = sin  0 sin  d> 
d3  = a.8  0 

XP4  = XP3  * b3 

YP4  = YP3  * 83 

ZP4  = ZP3  * h3 


Rear ranginp, 


cos  4>  = [(Xp4  tan  013  (d2C23  + d3C33)  - Yp4(d2C23  + d3C33)  (A-10) 

+ ZP4  ^d2c22  + d3c32)l/^zp4  tan  a3  sin  0 ) 

Equations  (A-9)  and  (A>10)  are  combined  and  the  angles  <f>  and  6 are 
determined  explicitly. 


Mi(Yp-b) 


zP4tana3  + 75yIT 


°i  - r^r 

P 


4)  = tan 


M tan  a,  + 

P 1 2 


(A- 1 1 ) 


0 = tan 


where 


XP~a 

(Xp-a)  tan  ctj  sin  4>  -( Yp-b)  tan  a ^ cos  4>  + Zp  cos 

Mj  = -C33Yp4  + C32Zp4  - C33Xp4  tana3 
M2  = *C23Yp4  + C22ZP4  + C23xP4  Un  a3 


T (A-  1 2) 


Thus,  Equations  (A-6),  (A-7),  (A-8),  A-ll),  and  (A- 12)  form  a complete 
set  of  formulae  for  determining  positions  and  angles  irom  Cameras  1 
and  3. 
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2. 


Determination  of  Dipole  Position  and  Attitude  (Cameras  1 and  2) 


The  determination  of  dipole  position  and  attitude  from  the  two 
horizontally  aligned  cameras  is  accomplished  in  a manner  similar  to  that 
developed  for  cameras  l and  3.  The  results  are  described  in  the  fol- 
lowing text.  From  camera  1,  the  dipole  position  and  attitude  are 
described  by 


= slope  of  image  to  the  Y axis  in  the  YZ  plane. 

Similarly,  from  camera  2 the  dipole  position  and  attitude  is  given  by 

X = e2 
Z = c2 

°2  = slope  of  image  to  the  X axis  in  the  XZ  plane. 

The  dipole  position  with  respect  to  the  referenced  coordinates  is 


a[a^2~ni  U2“b)] 

a2-(£2-b)  (n^-b) 

afar^-f^  (n1-b)] 
a2-(£;2-b)(n1-b} 


(A- 1 3) 


(A-  14) 


(,-b)  lU-blUt  « C?)  - (t;,^  ♦ V-,)) 
2 a2-K2-b)(n1-b) 


(A  - 1 5) 


The  dipole  orientation  is  obtained  from 


$ ■ tan 


n.-b 

a 

tan  ♦ 

Cj/a 

Vb 

a 

tan  a2  ♦ 

C2/« 
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cos  $ (tan  a0 


tan  0o  tan  $ - 


It  should  be  noted  that  in  the  data  reduction  utilizing  cameras  1 
and  2,  poor  accuracy  of  orientation  can  occur  when  the  chaff  dipole  is  in 
a nearly  horizontal  plane  (when  crossing  the  Z = 0 plane).  If  at  this  time 
a j and  are  equal  to  or  nearly  zero,  then  small  errors  in  reading 
either  or  a,  result  in  large  changes  in  $ and  0 . The  other 
orthogonal  set  of  cameras  1 and  3 were  Incorporated  into  the  test 
facility  to  circumvent  this  problem. 

3.  Calibration 


The  test  facility  was  fabricated  originatty  with  a wire  grid  system 
in  the  XZ  and  YZ  planes.  The  integrated  intensity  of  light  on  the  grid 
wires  during  the  particle  descent  was  so  large  that  the  grid  often 
obscured  the  position  of  the  chaff  filaments.  After  determining  that 
there  were  no  noticeable  optical  anomalies  in  the  lenses  (except  for 
parallax),  the  grid  wires  were  removed  in  favor  of  a single  pendulum 
wire  with  three  bead  reference  marks,  each  12  inches  apart.  The 
center  bead  located  the  zero  intercept  for  the  XYZ  coordinate  system. 
The  three  cameras  were  located  in  the  following  positions  with  respect 
to  the  pendulum  grid; 


Camera  1 


X = a = 96  in. 
Y = b = 16  in. 
Z = 0.0  in. 


Camera  2 


X * b * 16  In. 
Y = a = 96  in. 
Z 2 0.0  in. 


Camera  3 X = b = 16  in. 

Y = d = 34. 4 in.  . 

Z 2 e 2 72  in.  I 

0O  = 16. 6 deg.  I 

The  calibrated  pendulum  was  used  as  the  tare  reference  in  the 
following  planes: 

Camera  1,  plane  YZ  , tare  1 = 12  in. 

Camera  2,  plane  XZ  , tare  2-12  in. 

Camera  3,  plane  X3Z3,  tare  3 = 5.  003  in. 
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4.  Computer  Program 

All  chaff  dipole  test  drops  are  analyzed  in  detail  using  a computer 
program  designated  ASTRO.  This  program  is  written  in  FORTRAN  4 
and  is  operated  on  a CDC  6600  system.  The  program  statements  and 
operation  are  described  in  Reference  5. 


| 


[ 

I 


I 

! 

j 


^ Mihora,  D.  J.  , Astro  Research  Corp,  , "Photogrammetry  of  Minute 
Chaff  Elements  in  a Quiescent  Test  Chamber,"  Astro  Report  No. 
ARC-TM-  320,  October  3.  1974. 
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APPENDIX  B 


AERODYNAMIC  DATA  REDUCTION  PROGRAM 


Aerodynamic  data  reduction  is  accomplished  by  computer,  using 
as  input  dipole  time-position-attitude  histories  and  dipoie  physical  char- 
acteristics data.  The  time-position-attitude  histories  are  obtained  in 
card  format  from  the  ASTRO  data  reduction  program.  Appendix  A, 
which  converts  the  photographic  imagery  to  a motion  sequence.  The 
principal  aerodynamic  data  reduction  operations  are  described  below. 

1 . Data  Input  and  Smoothing 

Step  1 The  raw  data  produced  by  ASTRO  is  read-in  and  trans- 
formed to  the  inertial  coordinates  X,  Y,  Z and  the  Euler  angles  ^ and 
0 , as  shown  in  Figure  B-i.  A time  base  is  established  from  the  rotat- 
ing shutter  sampling  frequency. 

Step  2 The  X and  Y horizontal  coordinates  are  translated  to 
inertial  coordinates  XQ  and  Y0 , with  their  origin  at  approximately_the_ 
center  of  the  dipole  descent  spiral,  by  subtracting  constant  values  (X,  Y) 
from  each  point  of  the  run.  The  values  (X,  Y)  are  determined  from 
initial  screening  of  the  data. 

Step  3 Rectangular  position  coordinates  XG,  Y0  are  converted 
to  polar  coordinates  R and  V as  shown  in  Figure  B-2,  in  preparation  for 
determining  X0.  Y0,  X0,  and  Y0. 

v = tan*1  (Y0/X0)  (B-l) 

R = X0/cos  v (B-2) 

Step  4 The  variables  R,  v , Z , \p  and  0 are  locally  smoothed 
using  a quadratic  least- squares  routine  with  time  as  the  independent 
variable.  Optionally,  five  or  seven  values  are  taken  for  each  smoothing 
interval,  with  the  mid-point  being  the  point  desired.  For  each  smoothing 
operation  a quadratic  equation  which  best  "fits"  the  raw  data  is  obtained, 
and  by  differentiation  of  this  equation  the  first  and  second  time  deriva- 
tives are  obtained. 
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X,Y,Z  are  inertial  axes 

x,y,2  are  body  axes  assuming 

no  rotation  about  x (i.e. 
fixed  plane  axes, 
y remains  In  inertial 
X-Y  plane) 


2 


Angle  of  Attack  and  Flight  Path  Angles 


Step  5 Velocities  and  accelerations  with  respect  to  the 
Cartessian  coordinates  (X0,  Y0,  XQ,  and  Y0)  are  obtained  from  the 
smoothed  values  of  R and  V and  their  derivatives,  as  follows: 


= ft  cos  v - R v sin  v (B-3) 

Y0  = ft  sin  v - R v C03  v (B-4) 

£0  * 5 coe  v -2ft  v sin  V -R  v2  cos  v -R  v sin  v (B-5) 

VQ  = R sin  v *2ft  v sin  v -R  v2  sin  V +R  v cos  v (B-6) 


Step  6 The  inertial  velocities  and  accelerations  are  trans- 
formed to  fixed-plane  velocities  (u,  v,  w)  and  accelerations  (u,  v,  w) 
using  the  usual  Euler  transformation  matrix  (with  $ * 0 ). 


u,  A 

cos  ^ cos  0 sin  cos  0 -sin  6 

r vv 

V.  V 

■ 

- sin  tp  cos  0 0 

*o.  *o 

.V,  *. 

cos  i>  sin  6 sin  sin  6 cos  6 __ 

. io.  5o  . 

Step  7 The  total  angle  of  attack,  a , the  angle  of  attack  plane 
orientation  angle,  £ , the  flight  heading  angle,  T , the  desce  it  flight  path 
angle,  y , the  orientation  angle,  J , and  the  total  velocity.  Vj  . are 
obtained  according  to  the  following  equations.  The  angles  are  described 
in  Figures  B-3  to  B-5, 


a ■ tan-*  [Vv*  ♦ v2/u]  (B-b, 

£ - tan-1  (w/vl  (B-9) 

Y - tan-1  IY0/X0)  (B-10) 

y - tan-1  [ z/Vxq  ♦ *o  ] (B-ll) 

4>  ■ tan-*  £(u*v)/(v»Vj)J  ( B - 1 2 ) 

VT  - V u2  ♦ v2  ♦ v2  (B-13) 
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3. 


Aerodynamic  Forces 


Step  8 The  aerodynamic  forces  in  the  fixed-plane  coordinate 
system  are  determined  from  the  following  equations: 


Fx  = m (u  + g sin  9 ) 
Fy  = m v 

Fz  - m (w  - g coi  9 ) 
where  in  is  *he  dipole  mass  nnd  u, 


(B-14) 
( B- 1 5) 
(B-16) 

w are  obtained  from  step  ((>). 


Step  9 Using  the  angle  C , the  Fy  and  Fz  forces  arc  resolved 
into  force  components  in  the  pitch  plane  (angle  of  attack  plane)  denoted  by 
Fjq  or  normal  force  and  perpendicular  to  the  pitch  plane,  denoted  by  Fg 
or  side  force.  The  forces  Fjq  and  Fg  are  t'n»*n  converted  to  coefficient 
form  in  accordance  with  the  following  equations.  The  transformation  is 
described  in  Figure  B-6. 


cN 


CSF 


1?  ’ 5s  ( ' Fy  c°'  c ' F*  *'"  5) 
5?  ‘ 5s  ( Fy  ,in  E • F»  c"*  c ) 


( B- ) 7) 


( B-  1 8) 


where  q is  the  dynamic  pressure,  1/2  p V^,  , and  S is  the  aerodynamic 
reference  area,  taken  to  bo  the  dipole  maximum  projected  area,  i.  c.  , 
length  times  width.  The  axial  force  coefficient  is  obtained  similarly 


CA 


(B-19) 


4.  Aerodynamic  Moments 

Step  10  The  fixed-plane  moments  arc  given  by  the  following 
equations,  with  the  intrinsic  assumption  that  the  dipole  transverse 
moments  of  inertia,  ly  and  Iz  are  equal  (ly  = Iz  = I). 

Mx  = I p 
x xr 


( B-20) 


Figure  B-7.  Moment  Vector  Definitions 
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The  dynamic  pressure,  q , and  the  aerodynamic  reference  area,  S , are 
aa  before;  l is  the  aerodynamic  reference  length,  taken  to  be  the  dipole 
length. 


VI 
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APPENDIX  C 

CHAFF  AERODYNAMIC  THEORY 


1.  Theoretical  Aerodynamics  of  Symmetric  Dipole  at  Low  Reynolds 

Number 

The  complete  motion  of  bodies  with  poor  aerodynamic  shape  can- 
not, v/ithin  the  present  state  of  fluid  mechanics,  be  calculated  from  the 
equations  of  motion  for  the  flow.  Therefore,  a theoretical  treatment  of 
chaff  aerodynamic  behavior  must  be  restricted  to  those  fluid-dynamic 
conditions  which  are  amenable  to  analysis.  Thus,  we  find  that  the  aero- 
dynamic theory  is  limited  to  flows  for  which  the  viBcous  effects  can  be 
accounted  for  in  a suitably  approximate  manner,  and  to  body  shapes  and 
geometries  which  are  consistent  with  allowable  boundary  conditions. 


The  flow  past  a small  body,  such  as  a chaff  element,  can  usually 
be  described  bv  parameters  U,  , U , P , and  p^,  which  are  the  free 
stream  velocity,  characteristic  length,  fluid  viscosity,  fluid  density,  and 
fluid  pressure,  respectively.  By  dimensional  analysis  the  number  of 
parameters  may  be  reduced  to  one,  namely,  the  Reynolds  number,  Rj«j  , 
defined  by 


or 


U ft 
v 


where  v is  the  kinematic  viscosity. 

Before  taking  up  specific  cases  it  is  useful  to  consider  the  basic 
characteristics  of  fluid  flow  with  increasing  Reynolds  number.  There 
are  three  principal  regimes,  which  are  more  or  less  distinct:  1)  steady 
laminar  flow,  2)  a regime  in  which  regular  time-dependent  fluctuations 
occur  in  the  flow,  and  3)  a regime  in  which  part  of  the  flow  is  turbulent. 
The  motion  of  chaff  involves  only  the  first  two  of  these  three  regimes. 

The  solution  of  the  exact  Navier-Stokes  fluid-dynamic  equations 
for  the  conditions  of  R jq  lending  to  zero  is  due  to  Stokes.  A character- 
istic property  of  these  solutions  is  that  the  inertial  forces  of  the  fluid  are 
neglected  compared  to  the  v.scous  and  pressure  forces.  For  three- 
dimensional  flow,  solutions  to  the  Stokes  equations  may  be  found  for 
specified  boundary  conditions  on  a solid  and  at  infinity.  The  Stokes 
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resistance  formula  for  a sphere  of  radius,  a , is  (see  for  example. 
Reference  6). 

D _2*f  t 2 Ua 

CD  = i/2pu2  na2  : «N  ’ RN  t T"  (C' 

Similarly,  for  a disk  one  obtains  the  case  of  the  flow  perpendicular  to 
the  stream 


D = 16  U Ua  (C-3) 

and  for  the  disk  edgewise  into  the  flow 

D = -y  u Ua  (C-4) 

These  results  emphasize  one  of  the  important  features  at  low  Reynolds 
number:  the  resistance  is  relatively  insensitive  to  body  orientation. 

Although  drag  predictions  using  Stoke' s equations  have  been  found 
to  be  accurate  at  sufficiently  low  Reynolds  numbers,  difference  between 
experimental  data  and  the  Stokes  predictions  become  significant  when  the 
Reynolds  number  attains  a value  near  unity.  Also,  because  of  the  inher- 
ent assumptions  in  the  Stokes  solution,  the  Stokes  theory  cannot  deter- 
mine the  resistance  of  two-dimensional  bodies,  such  as  would  be 
representative  of  long  chaff  elements. 

The  so-called  Oseen  equations  are  obtained  by  linearization  of 
the  Navier-Stokes  equations  about  zero  velocity,  and  take  into  account 
viscosity  transport.  The  Oseen  equations  are  an  improvement  over  the 
Stokes  equations  because  they  provide  for  approximations  which  are 
uniformly  valid  over  the  entire  flow  field  at  low  Reynolds  numbers,  and 
because  they  more  accurately  represent  the  physical  flow,  i.  e.  , a wake 
in  the  lee  of  the  body  is  predicted.  The  Oseen  equations  have  been  used 
as  a starting  point  for  the  analysis  of  the  flow  around  various  shapes,  as 
will  bo  discussed  in  the  following  sections  on  Drag  and  Lift.  Here  we 
shall  consider  only  the  flow  characteristics  of  two-dimensional  circular 
cylinders  of  radius,  a . The  most  comprehensive  analysis  of  this  case 


Lagerstrom,  P.  A.  , Laminal  Flow  Theory,  "Theory  of  Laminar 
Flows,  " Vol.  IV.  High  Speed  Aerodynamics  and  Jet  Propulsion, 
Princeton,  New  Jersey,  Princeton  University  Press,  1964. 
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la  due  to  S.  Tomotika  and  T.  Aoi.^  The  drag  coefficient  ia  found  to 
be  of  the  form 

CD  = ^ ( « + 0 (Rn  e )J  (C-5) 

where  e is  the  flow  perturbation  parameter. 


Tomotika  and  Aoi  show  that  by  expansion  of  the  stream  function 
the  first  order  prediction  for  the  cylinder  drag  coefficient  is 


CD  = 


iLE  * , Y = Euler' b constant 

RN  1/2  - Y - £n  (Rn/8) 


Rn  (2.  002  - £n  RN)  (C‘6i 

Tomotika  and  Aoi  also  show  that  the  drag  of  the  circular  cylinder,  as 
determined  from  Oseen's  equations,  is  composed  equally  of  pressure 
and  frictional  drag  components. 

To  appreciate  the  significance  of  the  resistance  formulus  for 
cylinders,  it  is  appropriate  to  examine  the  effect  of  cylinder  radius 
(i.  e.  , chaff  3ize)  on  both  Reynolds  number  and  the  drag  coefficient. 
Unfortunately,  the  relationship  between  Cd  and  the  cylinder  radius,  a , 
is  given  by  a transcendental  function  which  must  be  solved  by  iteration. 
Also,  to  accomplish  the  solution  it  is  necessary  to  consider  the  flight 
velocity,  which  for  the  cylinder  in  steady  vertical  descent,  is  given  by 


U 


Tia  w 
P CD 


(C-7 


where  w is  the  chaff  weight  per  unit  volume.  Combining  Equations 
(C-l),  (C-6),  and  (C-7)  we  obtain  the  desired  relationships  between  Cd, 
Rn*  and  a*  Figure  C-l  shows  the  resulting  var  iation  of  Rn  with 
cylinder  radius  for  values  of  a between  0,0001  and  0.1  inches,  with 
w = 0.075  - 0,10  lb/in^  and  p and  p evaluated  for  standard  sea  level 


7 

Tomotika,  S,  , and  T.  Aoi  (University  of  Kyoto,  Japan),  "The  Steady 
Flow  of  Vi scous  Fluid  Past  a Sphere  and  Circular  Cylinder  at  Small 
Reynolds  Numbers,  " Quart.  Journ.  Mech.  and  Applied  Math.  , 

Vcl.  Ill,  Pt.  2(1950). 
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conditions.  Inspection  shows  a cylinder  radius  of  about  1 x 10"^  inches 
(2  mil  dipole  diameter)  corresponds  to  a Reynolds  number  of  about  0.55 
for  w = 0.1  lb/in'.  The  drag  coefficient  for  this  configuration,  in 
accordance  with  Equation  (C-6),  is  Cq  = 17.6. 

Thus,  with  respect  to  the  principal  flow  regimes,  it  is  significant 
that  a Reynolds  number  of  unity  marks  not  only  the  upper  extreme  of  the 
steady  laminar  flow  regime  but  also  the  approximate  maximum  Reynolds 
number  which  will  be  experienced  by  a cylindrical  dipole  with  diameter 
of  0.  002  inches  or  less.  A Reynolds  number  of  unity  is  also  the  approx- 
imate limit  for  which  the  Oseen  solutions  are  valid. 

At  Rj\j  = 40  the  previously  mentioned  second  flow  regime  begins 
with  the  appearance  of  the  phenomenon  of  vortex  shedding.  Eddies, 
which  at  lower  values  of  Rjq  had  become  distended,  now  periodically 
breaks  away  from  the  cylinder  and  flow  downstream  in  the  wake.  The 
nondimensional  shedding  frequency  is  given  by  the  Strouhal  number,  S, 
which  is  defined  as 

Strouhal  number  = S = (C-8) 

where  n is  the  shedding  frequency,  and  d the  cylinder  diameter. 

The  regime  10  < Rjq  < 10*  is  also  discussed  by  Lager strom^). 
Foil  chaff  elements,  which  have  a greater  characteristic  dimension  than 
cylindrical  dipoles  (i.  e.  , strip  width)  have  typical  descent  Reynolds 
numbers  of  the  order  of  1.0  - 100.0,  and  thus  can  be  expectea  to  experi- 
ence flow  conditions  representative  of  the  second  flow  regime.  In  this 
range  of  Reynolds  number  experimental  results  described  byLagcrstrom 
show  that  the  total  drag  coefficient  of  the  cylinder  continues  to  decrease 

with  increasing  Rjq  . with  Cqj  approximately  equal  to  4/  >/ Rjsj . The 
ratio  of  pressure  drag  to  friction  drag  at  Rj«q  = 40  increases  to  about 
1.  5,  whereas  at  Rjq  <<  1 the  ratio  is  unity.  While  the  above  comments 
apply  primarily  to  cylinders,  similar  characteristics  are  noted  for  bluff 
shapes  such  as  flat  plate,  etc.  For  the  flat  plate  at  zero  incidence  the 
viscous  drag  can  be  evaluated  from  boundary  layer  theory.  Using  the 

Blasiua  profile  the  overall  plate  drag  is  Cf  = 1.328  y/  Rjq  , which  is  in 
good  agreement  with  experiments  in  the  Reynolds  number  range  under 
discussion. 

The  foregoing  discussion  has  been  concerned  primarily  with  fluid 
flow  and  the  body  drag  force.  The  lift  force  and  moment  on  an  inclined 
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axlsymmetric  body  in  low  Reynold*  number  flow  is  of  considerable  inter- 
est with  respect  to  chaff  dynamics.  Unfortunately,  solutions  for  the  lift 
and  moment  have  not  received  as  much  attention.  Cox**^  has  found  that 
the  Cseen  equations  do  not  in  general  furnish  an  exactly  correct  value  of 
the  lift  force.  By  the  use  of  the  singular  perturbation  method  more 
accurate  functional  representations  of  the  aerodynamic  forces  to  OlRfj) 
and  0(R^j  In  Rjsj)  are  achieved.  However,  useful  solutions  for  the  lift 
force  on  two  dimensional  elliptic  cylinders  and  plates  with  arbitrary 
incidence  have  been  obtained  from  the  Oseen  equations  bv  Tomotika  and 
Aoi'7',  accurate  to  O(Rjsj).  From  their  results  it  is  found  that  the  drag 
increases  slightly  with  an  increase  of  either  thickness  or  angle  of 
incidence,  and  that  the  lift  decreases  with  thickness  ratio  while,  as  a 
function  of  incidence,  it  has  a maximum  at  about  45  degrees. 

The  aerodynamic  moment  on  an  inclined  body  at  low  Reynolds 
numbei  it  at  most  small.  Willmarth  et  al,  in  their  investigation 

of  freely  falling  disks  refers  to  the  work  of  Gans,  who  was  able  to  show 
that  for  low  Reynolds  number  creeping  motion  there  is  no  aerodynamic 
moment  produced  by  fluid  pressure  on  a body  that  has  three  mutually 
perpendicular  planes  of  symmetry.  For  higher  Reynolds  numbers  no 
direct  computation  of  the  aerodynamic  moment  as  a function  of  body 
incidence  has  been  reported;  thus,  one  is  ted  to  consider  more  general 
analyses  of  body  angular  motion. 

The  angular  motion  of  a solid  body  in  a fluid  can  be  treated  from 
various  theoretical  viewpoints.  The  motion  of  a body  in  an  ideal  fluid  is 

g 

Cox,  R.  G.  , "The  Steady  Motion  of  a Particle  of  Arbitrary  Shape  at 
Smalt  Reynolds  Number,"  Journal  of  Fluid  Mechanics  (1965),  Vol.  23, 
Part  4,  pp  625-b43. 

9 _ 

Tomotika,  S.  and  T.  Aoi  (Univ.  cf  Kyoto  Japan),  "The  Steady  Flow  of 
a Viscous  Fluid  Past  an  Elliptic  Cylinder  and  a Flat  Plate  at  Small 
Reynolds  Numbers,  " Quart.  Journ.  Mech.  and  Applied  Math.  , Vol. 

VI.  Pt.  3 (1953). 

10 

Will  narth,  W.  E.  , et  al,  "Investigations  of  the  Steady  and  Unsteady 
Motion  of  Freely  Falling  Disks,"  Aerospace  Research  Laboratories 
Report  No.  ARL  63-176,  October  1963,  AD  420  913. 
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discussed  at  length  by  Lamb^  and  reference  la  made  to  the  earlier 
work  of  Kirchhoff,  who  showed  that  in  the  case  of  the  ellipsoid,  stable 
motion  is  possible  only  when  the  direction  of  motion  la  along  the  least 
axis.  For  viscous  fluids  Cox  obtains  a similar  result,  showing  that  the 
motion  of  the  oblate  spheroid  (discoid  object)  is  stable  only  when  the  flow 
is  parallel  to  the  axis  of  rotational  symmetry  and  that  the  prolate 
spheroid  (needlelike  object)  is  stable  or>’  r when  the  flow  is  perpendicular 
to  the  axis  cf  rotational  symmetry.  For  creeping  flow  (i.  e.  , very  low 
Reynolds  number)  Happel  and  Brenner^^  have  applied  symmetry  con- 
siderations to  the  equations  for  the  hydrodynamic  stress  tensor,  and 
have  determined  the  motion  characteristics  of  various  classes  of  bodies. 
The  results  show  that  ellipsoids,  or  bodies  of  revolution  with  fore-aft 
symmetry,  which  are  of  uniform  density,  can  attain  rotationally  free 
terminal  descent  for  all  possible  orientations.  This  has  bem  verified 
by  experimental  observation  of  small  spherically  isotropic  bodies,  such 
as  tegular  polyhedrons,  which  are  noted  to  fall  without  rotation. 

Happel  and  Brenner  also  discuss  the  motion  (and  force  tensor)  of 
needlelike  bodies.  The  results  are  of  interest  because  of  the  applica- 
bility to  chaff.  If  the  body  is  a needle-shaped  ellipsoid,  with  the  (3) 
axis  denoted  as  the  axis  of  rotational  symmetry,  then  the  translation 
resistance  tensor  is 


Ki  = K2 


‘ 8 v <>  *] 

, In  (2  <p)  ♦ 1/2  J C ’ 


*n  (2  $)  - 1/2 


PKU 


(C-9) 


where  $ = £/2c  and  l is  the  length  of  the  needle  and  c its  radius. 

For  large  t/c,  (where  the  constants  in  Equation  (C-9)  can  be  neglected 
in  comparison  with  the  logarithmic  terms)  the  flight  path  of  the  needle, 

* * Lamb,  Sir  Horrace,  "Hydrodynamics,"  Sixth  Edition,  Dover 
Publications,  New  York, 

12 

Happel,  John,  and  Howard  Brenner,  "Low  Reynolds  Number 
Hydrodynamics,"  Prentice-Hall  International  Series  in  the  Physical 
and  Chemical  Engineering  Sciences,  Prentice-Hall,  Inc., 

Englewood  Cliffs,  N.  J. 
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descending  under  the  influence  of  gravity,  has  been  computed.  The 
angle,  y , between  the  vertical  and  downward  direction  of  motion  is 


Y 


- gjn  2 6 
3 + cos  2 6 


} 


(C-10) 


where  8 is  the  orientation  of  the  needle  in  accordance  with  the  following 
sketch. 


The  glide  angle,  y , is  maximum  when 


8 


54.8  degrees 


(C-ll) 


and  corresponding  direction  of  movement  is  given  by 

V = tan  * (0.354)  = 19.5  degrees  (C-12) 

A more  general  result  can  be  obtained  for  the  glide  angle,  y , 
using  the  ratio  of  the  transverse  and  axial  resistance  factors.  Thus,  if 
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we  define  an  axial  resistance  factor 


K3  = 


and  transverse  resistance  factor 


and  their  rati  > as 


we  obtain 


‘1.2  “ MU, 


*1,2 


y = tan 


■ 1 I - sin  2 0 


which  is  maximum  for  0 


•5--' 1(^)1 


(C-13) 


Aerodynamic  Forces  Acting  Upon  a Symmetric  Dipole 


Cylindrical  Elements  - Flow  Perpendicular  to  Axis  For  cylin- 
drical dipole  with  radii  less  than  0,001  inches,  the  Reynolds  number  in 
steady  descent  will  be  If  ss  than  unity,  and  the  drag  coefficient  as  deter- 
mined from  solution  of  tne  Oseen  equations,  l.  e.  , Equation  (C-6),  is 
applicable  when  the  chaff  element  has  near  horizontal  orientation.  Since 
this  result  is  two  dimensional,  a correlation  should  be  applied  for  finite 
length  cylinders.  Although  no  theory  is  available  for  short  cylinders, 
we  can  use  the  results  of  Burger  as  presented  by  Happel  and  Brenner 
for  estimating  the  effect  of  the  slenderness  ratio.  Burger's  drag  rela- 
tion for  a cylinder  moving  perpendicular  to  its  axis  is 


F = - 


8 it  p 


Hn  (2a/b)  + 0.5 


(C- 14) 


where  2a  is  the  cylinder  length  and  b the  radius.  The  drag  coefficient 

2b  U 

based  on  area  4ab  and  Rjq  = — - — is 


fc. ' A .1  * ; ^ 'fl * '.yS.  * W **» ■$,>*«*, 


l( ! \ 

N \ tn  ( 2a/b)  + 0.5  ) 


(C-15) 


As  an  example,  the  effect  of  chaff  length  on  drag  ia  determined  uaing 
Equation  (C-15).  Consider  the  ratio  of  drag  coefficient  for  1-mil 
diameter  chaff  of  lengths  0.25  and  2.00  inches.  From  Equation  (C-15) 
we  obtain  the  result 


CD  (.25"  length) 
CD  (2.00"  length) 


= 1.35 


Cylindrical  Elements  - Flow  Parallel  to  Axis  The  determina- 
tion of  the  drag  of  cylindrical  chaff  elements  in  axial  flow  (or  at 
incidence)  requires  methods  which  differ  from  those  used  for  calculation 
of  the  cross  flow. 

In  the  case  of  needlelike  bodies  in  axial  flow,  the  boundary  layer 
near  the  rear  of  the  body  is  large  compared  to  the  body  radius,  while 
near  the  nose  of  the  boundary  layer  thickness  is  small.  Glauert  & 
Lighthlll^ 1 i'  suggest  a mixed  solution  which  provides  valid  values  of 
the  skin  friction  for  all  values  of  the  parameter  Vx/Ua^,  where  x is 
the  distance  from  the  front  of  the  cylinder.  The  total  momentum  defect 
(which  is  proportional  to  the  drag  coefficient)  is  presented  in  graphical 
form  in  Figure  C-2. 

A more  approximate  solution  for  the  axial  drag,  due  to  Burger, 
is  presented  by  Brenner; 


4 n u Ua 


In  ( 2a/b)  - 0.72 


(C-16) 


where  the  cylinder  length  is  2a  and  the  radius  is  h.  The  drag  coeffi- 
cient, based  on  area  4ab,  and  R^  = U2a/v  is 


4 n / a/b \ 

RN  \ *n  (2a/b)  - 0.72  / 


(C-  17) 


Glauert,  M.  B.  and  M.  J.  Lighthill  F.  R.  S.  , "The  Axlsymmctric 
Boundary  Layer  on  a Long  Thin  Cylinder,  " Proceedings  of  Royal 
Society  of  London,  Series  A,  Mathematical  and  Physical  Sciences, 
12  July  1955. 
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cylinder  radius 
cylinder  length 


Figure  C-2.  Drag  Parameters  of  Long  Cylinder  in  Axial  Flow  (Theory  of  Glauert  and  llghthlll) 


K-  ■ \ asm* 


To  compare  Glauert’s  and  Burger'*  results,  drag  coefficients  have  been 
computed  for  a 1 -mil-diameter  dipole  2-inches  in  length,  descending  in 
vertical  orientation  at  U = 1 ft/sec  at  sea  level.  We  obtain 


Glauert's  method 
Burger's  method 


CD  - 7.  50 
CD  = 3.  14 


Apparently,  at  U - 1 ft/sec  Burger's  approximate  solution  based  on 
creeping  flow  is  unreliable. 


Strip  Elements  - Flow  Perpendicular  to  Axis  The  drag  of  a 
two-dimensional  plate  (or  elliptic  cylinder)  with  axis  perpendicular  to 
the  flow  has  been  determined  by  Imai(^\  using  a general  solution  of 
Oseen's  equation  involving  successive  approximations  in  powers  of  the 
Reynolds  number.  The  drag  of  the  plate  at  zero  incidence  is 


white  for  90  degrees  incidence 


(C- 18) 


\ 


14 

Imai,  Isao,  (Univ.  of  Tokyo,  Japan),  "A  New  Method  of  Solving 
Oseen's  Equations  and  Its  Application  to  the  Flow  Past  an  Inclined 
Elliptic  Cylinder,"  Proceedings  of  Royal  Society,  A 224  (1954), 
pg.  141. 


J 03 


where 


(C- 19) 


and  Rjvj  is  based  on  the  width  of  the  plate,  2a  , 

Lift  (Normal  Force)  For  strip  chaff  there  are  two  sources  of 
aerodynamic  lift;  1)  the  inclination  of  the  longitudinal  axis,  and  2)  the 
orientation  of  the  plate  about  its  axis.  These  are  distinguished  in  Figure 
C-3. 


A theoretical  solution  for  the  effect  of  a on  the  lift  and  drag  of 
the  strip  type  element  is  not  presently  available.  For  a = tt  /2 , i.  e.  , 
with  pure  cross  flow,  Imai  gives  both  Cl  and  Cq  tor  the  two-dimen- 
sional flat  plate  as  a function  of  <t>.  The  resulting  variations  of  Cl 
and  Cq  with  4>  for  three  representative  Reynolds  numbers  are  shown 
in  Figure  C-4.  The  variation  of  Cl  with  $ is  teen  to  be  closely 
approximated  by  Cl  = (^l^MAX  ®'n  ^ * • 

For  a t v/2  cross  flow  theory  can  be  invoked.  The  normal 
force  can  be  computed  from  the  cross  flow  component,  Ujsj  , and  the 
cross  flow  Reynolds  number.  The  axial  force  can  be  computed  from 
either  the  akin  friction,  or  from  the  Stokes  solution  for  the  equivalent 
cylinder. 

As  a first  approximation  to  the  cross  flow,  we  can  take  the 
leading  terms  in  Imai's  formulas  and  obtain  the  total  normal  force 
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coefficient  as 


Cfj  ( a,  $) 


where  S 


(C-20) 


rN  * 


c = width  of  strip 
U = total  velocity. 


Thus,  the  normal  force  consists  of  one  component  which  is  independe  it 
of  ♦ , and  a second  component  which  depends  upon  both  a and  4>  . 

Side  Force  In  the  case  of  a foil  dipole  there  exists  a side  fi  rce 
as  well  as  a normal  force  and  axial  force  (see  Figure  C-3).  The  sid2 
force  coefficient,  Cgp  . is  derived  in  a manner  similar  to  that  used  for 
the  normal  force,  with  the  result 


cSF(3.*>  ■ (^)(-=f) 


sin  a • sin  2 # 


(C-21) 


This  force  is  seen  to  be  identical  in  magnitude  to  the  ♦ dependent 
normal  force,  but  with  tt/2  difference  in  phase. 

3.  Aerodynamic  Characteristics  of  an  Asymmetric  Dipole 

Dipole  asymmetries  provide  the  driving  force  for  all  spiral  and 
glide-type  motions.  Asymmetries  may  be  either  inertial  or  geometric 
in  form.  Several  types  of  asymmetry  will  usually  exist  simultaneously, 
however,  the  magnitude  of  each  type  of  asymmetry  will  vary  widely 
between  different  dipoles. 

All  asymmetries  are  referenced  to  a basic  dipole  configuration 
which  is  assumed  to  have  at  least  three  planes  of  mirror  symmetry.  In 
addition,  it  is  assumed  that  all  glass  dipoles  have  basic  axial  symmetry. 
Representative  symmetric  dipoles  are  sketched  on  the  following  page. 


i 

\ 
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a) 


Glass-Type  Dipole 
(axially  symmetric 
about  x ) 


bl  Foil-Type  Dipole 
(mirror  symmetric) 


Center -of -Mass  Offset  The  most  elementary  chaff  asymmetry 

is  a displacement  of  the  mass  center  from  the  origin  of  the  axes  of  sym- 
metry (see  sketch  below). 


This  type  of  mass  asymmetry  results  in  aerodynamic  moments 
about  the  center-of-mass  of  the  form 

AL  = - Fz  Ay  + Fy  Az 

AM  * Fz  Ax  - Fx  Az  (C-22) 

AN  . ■ - Fz  Ax  + Fx  Ay 

where  Fx,  Fy,  Fz  are  the  aerodynamic  force  components,  which  are 
unaffected  by  mass  asymmetry. 

Longitudinal  Bending  Longitudinal  bending  is  a significant 
source  of  aerodynamic  moment  when  the  dipole  axis  is  closely  aligned 
with  the  direction  of  flight.  This  is  because  the  local  viscous  force 
vectors  do  not,  in  general,  pass  through  the  cent' r of  mass.  A dipole 
with  longitudinal  bend  of  constant  bend  radius,  R , can  be  depicted 
rchematically  as  in  Figure  C-5).  The  moment  can  be  expressed  in  the 
form 


<IM  • F £ (6  - R)  cos  6 + R J ds 


(C-23) 


The  variable  is  the  viscous  force  per  unit  length  and  can  be  evalu- 

ated as  a function  of  s by  the  methods  of  Reference  13).  For  prelimin- 
ary analysis  the  average  viscous  force  can  be  used. 


' viscous 


2a 


(C-24) 


Substituting  Equation  (C-24)  into  Equation  (C-23),  letting  ds  = R d0  , 
and  integrating,  we  obtain 

„ . Fv„  [i  ,(i^B)sln(|j]  (C-25) 

or  in  coefficient  form 

°M  ■ C0.xlal(ll)  [!  <C-26> 

If  the  dipole  has  uniform  mass  density  and  if  R >>  a the  center-of- 
mass  offset  is  given  by 
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Thu* 


CM  - C°.*Ul  (fa)  [‘  + ( TS  ' 7 ) ■ln  f ] <C-28' 


A second  source  of  aerodynamic  moment  due  to  longitudinal  bend 
occurs  at  finite  angles  of  attack  and  is  done  to  the  variation  in  normal 
force  along  the  element.  (See  Figure  C-6). 

The  elemental  moment  due  to  can  be  expressed  as 


d CM  * 


2a 


(s) 


ds 


(C-29) 


Assuming  that  [Cjvj  (a)]  is  of  the  form 


N 


[CN(aj]  = a . (C-30) 

noting  that  a (s)  = aR  - 8 , and  letting  s = R 0 we  can  obtain 


M 


* TVJ'f  [*,tn(aR-  R )] 


ds 


( C - 3 1 ) 


C 

2V 


£sin  (a/R)  - (a/R)  cos  (a/R) J cos  nR 


The  Cj^  moment  due  to  the  normal  force  variation  (Equation  (C- 31 ))  is 
seen  to  oppose  that  due  to  the  axial  viscous  force  (Equation  (C-28)). 

Transverse  Bend,  Twist,  and  Kink  Examples  of  transverse 
bend,  twist,  and  kink  are  depicted  in  Figure  C-7. 

Since  these  shapes  are  not  amenable  to  direct  analysis,  the  func- 
tional relationship  between  the  aerodynamic  forces  and  moments  and  the 
dipole  attitude  and  configuration  geometry  must  be  established  largely 
from  symmetry  considerations.  To  this  end  it  is  useful  to  establish  for 
the  symmetric  dipole  force  ^nd  moment  coefficients  the  dependence  upon 
the  aerodynamic  variables  a and  t (see  Figure  C-3),  This  dependency 
is  noted  in  the  upper  half  of  Table  C-l.  For  the  symmetric  dipole  the 
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force  coefficient  relationships  for  CN  and  Cgp  are  consistent  with 
cross-flow  theory  (Equations  (C-20)  and  (C-21)),  while  the  moment 
coefficient  dependence  upon  a is  assumed  apriori.  The  harmonic 
coefficients  for  both  the  symmetric  and  unsymmetric  dipoles  are 
assigned  arbitrary  identities;  k^  denoting  the  force  coefficients  and 
Kj  denoting  the  moment  coefficients. 

For  asymmetric  dipoles  the  functional  dependence  of  the  forces 
and  moments  on  a and  4*  is  shown  in  the  lower  portion  of  Table  C-l, 
where  the  completeness  longitudinal  as  well  as  transverse  bend  is 
included.  The  formulas  for  transverse  bend  and  kink  are  postulated 
from  both  cross-flow  considerations  and  the  nature  of  the  configura- 
tional asymmetry.  In  the  case  of  twist  with  small  deflection  C , the 
functional  dependence  is  derived  analytically  from  the  force  and  moment 
definitions  for  the  symmetric  element. 

The  Kj's  in  Table  C-l  are,  in  general,  either  positive  or  nega- 
tive: thus,  if  one  or  more  moments  due  to  asymmetry  exist  for  each 
degree  of  freedom,  it  is  possible  for  a and  ♦ to  have  trim  values 
which  do  not  correspond  to  those  for  the  symmetric  dipole.  Also,  the 
moment  due  to  asymmetry  can  result  in  steady-state  angular  motions, 
wherein  the  moment  balance  is  completed  by  the  inertial  and  damping 
moment  contributions. 


4.  Aerodynamic  Damping 

Estimated  values  for  the  aerodynamic  damping  derivatives  corre- 
sponding to  the  plane  of  total  angle  of  attack  and  plane  of  the  side  moment 
are  required  for  analysis  and  prediction  of  dipole  angular  motion.  The 
damping  moment  contributions  are  particularly  significant  in  the  caee  of 
chaff  configurations,  not  only  because  the  damping  moments  tend  to  be 
large  as  a result  of  the  viscous  forces,  but  also  because  the  correspond- 
ing static  moments  tend  to  be  quite  small. 

Because  of  the  large  length  to  width  ratio  of  the  characteristic 
dipole,  it  is  appropriate  to  use  strip  theory  for  determination  of  the  local 
cross  force  along  a dipole.  It  is  assumed  that  the  local  force  has  the 
same  dependency  upon  flow  orientation  and  velocity  as  does  the  complete 
clement. 

Angle  cf  Attack  Plane  The  following  sketch  depicts  the  local  flow 
and  force  on  a dipole  undergoing  rotation  in  the  angle  of  attack  plane. 


/ 
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It  is  assumed  that  the  local  normal  force  coefficient  for  a nominal 
range  of  Reynolds  numbers  can  be  expressed  as 


CN  (a  , V) 


sin  a 


(C- 


where  is  a constant  corresponding  to  the  normal  force  on  a com 

plete  dipole. 


The  local  damping  moment,  M(q),  is  now  given  by 


dM 


(vfsln  “)  JP  Vl 


c x dx 


(C-; 


vi  ■ 


where 


V - qx  sin  a 


I 


1 


The  damping  derivative  has  the  same  inverse  dependence  upon  velocity 
as  do  the  force  coefficients. 

Side  Moment  Plane  For  the  side  moment  plane  the  local  cross 
force  and  velocity  components  are  depicted  on  the  following  page. 
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The  local  damping  moment  N(r)  is  given  by  the  differential 


dll  * [CN  (a,  V)  sin  O 2 p c x dx 

where 

. r rx 

sin  C = TT- 

vn 

Vfj  * V sin  a 

Vi  * >J  V2  + (rx)2 

Again  using  the  relation 


and  letting  Vj  = V (which  is  a reasonable  approximation  if  r is  small 
and  o does  not  approach  zero)  the  differential  moment  can  be  written 


dN 


— O 

Cjj  p c r r dx 
2 


(C-37) 


After  integrating  and  defining  the  damping  derivative  for  the  side  moment 
plane,  Cnr , as 


9 Cn 


where 

C n -*■ 

n 1/2  p c 


the  following  solution  for  the  damping  derivative  is  finally  obtained 


( C - 38 ) 


5,  Spiral  Dynamics 

The  characteristic  spiral  mode  of  a chaff  cipolc  results  from 

both  an  unsymmetric  trim  angle  of  attack,  , (i.  e.  , a^,  i 0,  it/2,...) 

and  an  unbalanced  static  side  moment,  SM,  which  causes  continuous 

» 

azimuthal  rotation,  <p  . A heuristic  model  of  the  aerodynamic  forces  is 
depicted  in  Figure  C-8  along  with  the  pertinent  motion  parameters.  The 
model  is  appropriate  for  those  cases  where  the  aerodynamic  roll  angle, 

♦ , is  either  constant  or  the  dipole  configuration  is  a;dsymmetric  and 
has  no  axial  spin. 

The  centrifugal  force  due  to  chaff  spiraling  motion  is  quite 
small*,  and  is  balanced  by  the  slight  inclination  of  the  normal  force 


* This  deduction  is  made  from  observation  of  dipole  spiral  motions. 
Comparing  the  calculated  centrifugal  force  (obtained  from  CF  = mR^| 
with  the  dipole  weight,  which  is  of  the  same  order  as  Fjq,  it  is  found 
*hat  CF  is  only  2-3  percent  of  Fjq  . 
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vector,  Fj^,  from  the  vertical  reference  plane  as  shown  in  view  A-A. 
The  inclination  of  Fjsj  is  possible  because  of  the  small  angle  of  yaw,  . 

For  steady- state  spiral  the  attitude  of  the  dipole  relative  to  the 
instantaneous  flight  path  remains  constant.  Thus  the  dipole  rotation 
about  itB  center-of-mass,  \ , is  also  constant  and  identical  to  the 
spiral  rate,  fl  . The  yaw  rate  can  be  computed  from 

.?  ^ Side  Moment 

w cos  0 » 

Yav  Damping  Moment 

which  in  coefficient  form  becomes 

1 CSM  /gV) 

* cos  0 = r " c^;  [rj 

where  and  Cnr  are  functions  of  Op  and  possibly  The 

spiraling  rate  ft  is 


Dispersion  Due  to  Steady-State  Spiral  The  radius  of  spiral,  R, 
is  most  readily  determined  from  the  kinematic  relationship 


(C-39) 


(C-40) 


R « 2-SSS-l 


where 


V 


is  the  total  velocity, 


/ 


2 W sin 
P Cp  S 


Y = 


is  the  flight  path  angle  and 


= lift  coefficient 
Cp  = drag  coefficient. 


(C-42) 


(C-43) 


(C-44) 


The  lift  and  drag  coefficients  are  functions  of  a and  possibly  ♦ . 
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The  horizontal  dispersion  of  a dipole  from  its  release  point,  r , 
depends  upon  the  descent  distance.  This  is  easily  seen  from  Figure  C-9. 
The  general  relationship  for  r is  given  by 


r = 2R  sin 


f — M 

\ 2R  CD  ) 


(C-45) 


The  maximum  radial  dispersion  is  readily  seen  to  be  r 
descent  distance  for  one  complete  spiral  is  obviously 


= 2R.  The 


, _ 2 R 

1 cr,/cn 


(C-46) 


For  h » hj  the  probable  value  of  radial  dispersion,  r , is  found  to  be 


r - -R 


(C-47) 
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APPENDIX  D 
ISOTROPIC  TURBULENCE 


This  appendix  presents  a method  of  digital  computer  modeling 
for  isotropic  three-dimensional  turbulence  applicable  to  chaff  dipoles. 
The  present  turbulence  model  has  greater  similarity  to  the  relation- 
ships for  heavy  particle  diffusion  than  to  the  "frozen-flow"  spectral 
models  common  in  aircraft  dynamics.  Existing  techniques  for  digital 
simulation  of  turbulence  (primarily  used  in  the  aircraft  industry)  are 
inappropriate  to  the  present  study  because  of  the  extremely  slow  descent 
velocities  of  chaff  dipoles. 

The  proposed  turbulence  model  for  chaff  dynamics  incorporates 
both  spacial  and  temporal  variations  of  the  turbulent  wind  vector,  and 
utilizes  as  input  actual  statistical  spectral  atmospheric  turbulence 
measurements.  A summary  of  the  dimensional  properties  of  atmo- 
spheric turbulence  parameters  is  given  in  Appendix  E. 

The  analyses  which  leads  to  the  final  set  of  equations  which  are 
used  in  the  digital  computer  simulation  are  presented  in  this  appendix. 
The  statistical  variables  and  functions  are  defined  in  Sections  D-l  and 
D-2;  the  joint  Lagrangian-Eulerian  correlation  function  is  presented  in 
Section  D-3;  Section  D-4  discusses  the  development  of  the  equations  for 
numerical  digital  simulation  with  the  6-DOF  computer  program. 


1.  Eulerian  Coordinates 


Two  distinct  reference  coordinates,  Eulerian  and  Lagrangian, 
are  commonly  employed  for  identification  of  turbulence  properties.  The 
Eulerian*  reference  is  used  in  the  overwhelming  majority  of  theories 
and  experiments,  for  example,  the  development  of  the  Navier  Stokes 
equations,  and  experimental  measurements  in  wind  tunnels,  and  from 
meteorological  towers.  Experiments  in  the  Lagrangian  coordinates  are 
difficult  because  they  involve  measurements  with  respect  to  the  mean 
fluid  movement.'  However,  developments  in  particle  dispersion  and  air 


* The  adjective  "Eulerian"  is  used  whenever  correlations  between  two 
fixed  points  In  a fixed  frame  of  reference  are  considered. 
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pollution,  including  atmospheric  measurements  with  tetroons  (Refer- 
ences 15  through  18), have  recently  yielded  a variety  of  Lagrangian 
statistical  properties. 

With  slowly  falling  particles,  such  as  chaff,  the  incremental 
changes  in  both  time  and  position  affect  the  ever-changing  turbulent  flow 
field,  as  seen  by  the  chaff  particle.  Herein  lies  the  principal  difficulty; 
that  of  combining  or  transferring  between  coordinates.  To  surmount 
this  problem  an  approximate  joint  correlation  has  been  developed  and  is 
presented  later  using  statistical  data  for  both  Eulerian  and  Lagrangian 
coordinate  systems. 


The  fluid  velocities  of  the  atmosphere,  U,  V,  W,  in  the  Euler 
coordinates  X,  V,  Z,  as  depicted  in  Figure  D-l,  are: 


u = xw  + U' 

(D-la) 

0 

where  Xw  = mean  wind  velocity 

and  U'  = turbulent  velocity  component. 

To  simplify  the  eventual  relationships,  the  atmospheric  mean  wind  is 
assumed  to  be  only  in  the  X direction.  Thus,  the  other  atmospheric 
velocities  are: 

> 

n 

> 

(D-lb) 

II 

3 

(D-lc) 

Assuming,  for  example,  a stationary,  random,  turbulent  field 
for  W'(t)  over  a time  span  of  t seconds,  the  spacial  power  spectral 
density  can  be  written: 

* (K)  • t“-  4k‘"o  srruio  f <"V^>  * <D-J> 

J0 

where,  in  position  space,  K(cycle/ft)  is  the  wave  number.  Similarly, 
a correlation  function  can  be  determined: 

t 

R (p)  - f w (Z,  t)  W (Z  + P,T)  dx  (ft2/sec2) 

^ o 
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where  p is  the  separation  distance.  T!».  correlation  coefficient  will 
be  denoted: 


R (p)  - V (Z)  V (Z  ♦ p)  ( D-3) 

The  over-bar  refers  to  stochastic  averages.  The  mean  value  of 
turbulence  is: 


V " V R (°°)  * 0 
The  mean  square  value  is: 

o2  « R (0) 
v 

Similar  results  are  obtained  for  the  U'  and  V'  velocities. 

Together  with  the  level  of  turbulent  fluctuations  ow  (rms),  the 
most  commonly  referenced  parameter  in  atmospheric  measurements  is 
the  intensity: 


l - ov/xw 


A dimensional  integral  scale  is  also  defined: 


(ft) 


(D-4) 


(D-5) 


In  the  investigation  (Reference  19)  of  high- Reynolds-number 
shear  flow,  a modified  integral  time  constant  on  Eulerian  space  is 
determined  to  be: 


I 


\ 


i 

I 

I 


I 


I 


T 


E 


dXv 

dZ 


+ 0/Lz 


(D-6) 


In  practice,  statistical  functions  using  the  cpacial  correlation  p 
are  rarely  considered.  Rather,  autocorrelation  functions  are  employed, 
and  time  averages  replace  ensemble  averages.  "Taylor's  hypothesis" 
(see  Reference  20)  is  commonly  referenced  for  this  conversion.  In 
this  hypothesis,  the  turbulent  flew  is  assumed  to  be  "frozen"  and  trans- 
ported with  the  mean  flow.  The  turbulent  components  do  not  change  with 
time  (which  Implies  that  diffusion  is  not  present).  The  frozen-flow 
modeling  suggests  that  ensemble  averages,  the  sequence  of  events  at  a 
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fixed  point,  are  nearly  equivalent  to  the  movement  of  an  unchanging 
pattern  past  this  point,  Also,  past  this  point  the  wave  number  K is 
related  to  spacial  frequency  by  velocity  W (see  Reference  15)  with: 

u>  » 2it  KW 

The  time  separation  variable  is: 

T * p/W 

The  correlation  function  (symmetric  about  the  origin)  and  the  spectra 
become: 


R (t) 


7 


♦ (u>)  COS  UT  da) 


Af  ’ 


(w) 


*/ 


R (x)  cos  (ut)  dT 


Comparisons  of  spectra  and  correlation  functions  for  the  atmosphere 
are  described  in  detail  in  Appendix  E.  The  simplified  correlation 
function  used  here  is  of  the  form: 


R (P) 


(D-7) 


Although  the  utilization  of  the  frozen-flow  model  is  quite 
universal  in  the  dynamic  analysis  of  flight  vehicles,  it  is  shown  in 
Reference  21  that  this  model  is  limited  to  vehicles  which  have  an  aero- 
dynamic velocity,  ua , va , wa  in  excess  of  approximately  40  , where 
0 is  the  rms  turbulent  velocity.  Typically,  the  maximum  aerody- 
namic velocity  of  chaff  in  the  Stokes  flow  regime  is 


va  = Z a 0 (a) 

and  thus  the  criterion  for  frozen  flow  are  not  satisfied. 


(D-8) 
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2.  Lagrangian  Coordinates 

For  chaff,  the  Eulerlan  modeling  must  be  augmented  with 
Lagranfian  terma  which  uescribe  the  time  variability  of  the  atmospheric 
wind. 

In  Lagrangian  apace,  an  integral  time  constant  is: 


w 

7“ 


Reference  23  indicates  that  this  time  constant  appears  related  to  the 
Eulerian  time  constant,  t£  . Therefore,  the  integral  time  scales  can 
be  related  by  a simple  constant  or  function  0 , where: 


VTE 


(D-9) 


Values  of  0 from  0.  3 to  10  have  been  reported  in  the  literature,  with 
values  of  1 to  6 most  prevalent.  For  turbulence  intensities  where 
i > 0.  15,  Reference  15  proposes  the  formulation: 

0 - 0.5/i 

Data  correlation  for  low  intensities  of  turbulence  at  an  altitude  of  2500 
feet  were  given  in  Reference  16  and  correspond  to  the  following: 


0 - e 


■6.5  i + 2.2 


(D-10) 


Reference  26  discusses  turbulence  behavior  at  high  Reynolds  numbers 
and  intenoities,  and  arrives  at: 

0 » h.O 

Reference  23  discusses  simple  relations  for  similarity  at  small  micro- 
scales and  arrives  at: 

0 a 1.0 


For  generalized  modeling  with  moderate  to  low  turbulence  intensities, 
Equations  (D-6),  (D-9),  and  (D-10)  are  combined  to  dete rmine  the 
Lagrangian  time  constant: 
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-6.5  i + 2.2 

*I£ 

dZ 


The  Lagrangian  correlation  function  (considered  in  References  23  - 25 
and  similar  in  form  to  Equation  (D-7)  is: 


I?  (t)  = 02  e"t//TL 


(D-12) 


3.  Joint  Correlation  Functions 

A joint  Lagrangian-Eulerian  correlation  function  is  now  defined 
which  incorporates  Equations  (D-7)  and  (D-12): 

R (t,p)  = o2  exp  - | [a  (t/TL)N  + b (p/Lg)1*  ] 1/,N  } 

Values  of  N = 2 and  a,  b = 1 are  chosen  which  yield  "smooth" 
functions  when  t or  p are  nearly  equal  to  zero  (i.  e.  , small  micro- 
scales). Reference  16  discusses  a similar  form  of  the  joint  correlation 
function  incorporating  a time  scale  of  3 =1.  It  is  assumed  that  this 
joint  correlation  function  is  aligned  with  the  mean  wind.  This  orientation 
is  acceptable  when  particle  velocity  is  significantly  larger  than  turbulent 
fluctuation.  However,  when  the  chaff  particle  velocity  is  of  the  order  of 
0,  the  aerodynamic  velocity  of  the  particle  ua,  va,  and  wa  may  be 
equally  large;  thus  any  orientation  (X,  Y,  Z)  might  be  denoted  as  the 
axis  of  principal  correlation  (i.  e.  , the  particle  becomes  lost  in  the  eddy 
dynamics).  However,  in  a statistical  sense  this  i'  of  little  note  when 
the  direction  of  principal  correlation  is  changed  wirh  time,  because  the 
eventual  average  direction  corresponds  to  the  Z direction  (i.  e.  , 

X - X*  a 0.  Y = 0,  Z i 0). 

Therefore,  the  principal  correlation  function  becomes: 

R£  (t,  Z)  = o2  exp  -|(t/TL)2  + (tZ/LE)2|l/2  (D-13) 


Secondary  correlation  functions  in  the  X and  Y directions  are 
also  required.  These  are  obtained  with  the  simplifying  assumption  that 
above  the  earth  boundary  layer  the  turbulence  is  isotropic  (Reference  27, 
p.  28).  For  isotropic  conditions  the  variances  of  the  turbulent  com- 
ponents are  all  equal,  i.  e.  : 
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( D- 14) 


From  Reference  28  (VonKarman  ) the  primary  and  secondary  correla- 
tion functions 


R*  (P)  » W'(Z)  W'(Z  ♦ p) 
Rn  (p)  - v*(z)  V'(Z  + p) 


V (z) 


V*  (z) 


W (Z  + p) 

are  shown  to  have  the  weil-known  relationship 

, 3 R»  (P) 

Rn  (o)  - R*  (p)  + ip  ~— 


V»  (Z  + p) 


(D-15) 


The  secondary  correlation  functions,  Equations  (D-16a),  (D-l6b)  are 
obtained  from  Equations  (0-13)  and  (D-15).  Equation  (D-13)  is  re- 
written as  Equation  (D-’.6c)  which  completes  the  following  set  of  corre- 
lation functions  for  XYZ  space: 


R (t,X) 


1 - 


t Z2 


8LX»- 


'i_+  (_i_f ' 

t2  u l x; 

Ltl 


2 11/2 


exp 


-b-w] 


1/2, 


(D-l6a) 
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R (t ,Y) 


R (t,Z) 


(D-  16b) 


( D-  16c) 


For  isotropic  flow,  the  scale  lengths  are  related  by  the  following  (see 
Reference  28,  Von  Karman), 


L%  B * DLy 


(D-17) 


4,  Numerical  Modeling 

Turbulent  flow  is  a random  process,  but  unlike  Brownian  move- 
ment the  fluid  displays  interaction  and  a memory  between  cells.  When 
sampled  over  long  time  periods,  turbulence  processes  display  Gau ssian 
number  averages  with  the  variance  oP  . The  following  generalized 
model  is  proposed  which  combines  memory  and  randomness,  and  accom- 
modates the  cor. elation  function  (Equations  (D-l6a,  b.  c)  ). 


u*  = 

m 

pu;.i 

+ v 1 - 

P2  Gm 

(D-  18a) 

V*  r: 
m 

+ V i - 

( D- 1 8b) 

W»  a 
m 

rw;-i 

+ -v/l  - 

R2  J_ 
m 

(D- 1 8c) 

The  values  P,  Q,  R are  memory  functions,  and  Gm,  Hm,  Jm  are 
random  numbers  chosen  from  a normal  distribution  of  zero  mean  and 
standard  deviation  aQ. 
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Equations  (D-18a)  through  (D-18~)  are  in  a particularly  attractive 
form  for  digital  modeling.  The  subscripts  m and  m-1  respectively 
designate  the  stepwise  present  and  previous  values  of  the  variable. 

The  following  stochastic  averages  are  formulated: 


U’ 

U’  , 

o 

* 0" 

p 

m 

m-1 

o 

U' 

U:  o 

* o2 

P2 

(D- 19) 

a 

m-2 

c 

u7* 

TFT 

o 

=>  0“ 

p3 

m 

m-3 

0 

• 

• 

The  correlation  function  which  identifies  the  previous  sequence 
is  of  the  form: 


R (t,  p)  = 0?  PN 


( D-20) 


where  N *»  t/At 

Then  the  correlation  function  is  simply: 


(t,  p)  = a2  e-log  (?)(t/At) 


(D-21) 


where  t = real  time 

and  At  = time  increment  used  in  6-DOF  computer  program. 

Memory  terms  P,  Q,  R can  now  be  determined  by  equating 
Equation  (D-21)  to  Equation  (D-16a),  etc.,  resulting  in  the  following 
expressions : 


(-“[-*]! 


( D-22a) 
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where 


■A 


where 


( D-22b) 


(D-22c) 


Equations  (D- 18a,  b,  c)  and  (D-22a,  b,  c)  were  programmed  into 
the  6-DOF  analysis  as  a subroutine  called  TURB1.  The  random  wind 
components,  G,  H,  J are  drawn  from  a subroutine  called  GAUSS. 
Currently,  the  values  of  Lx»  Ly.  an<*  > ant*  °X»  °Y»  anc*  °Z  are 
defined  for  the  isotropic  state,  and  are  given  by  Equations  (D-17)  and 
(D-14),  respectively.  The  program  initiates  with  the  three  turbulent 
wind  components  all  equal  to  zero. 

It  should  be  noted  that  the  equations  for  turbulence,  although  in- 
volving the  time  increment  At  and  chaff  velocity  Z , are  not  integrated 
directly,  but  constitute  a set  of  recursion  formulas  from  which  the  chaff 
aerodynamic  velocity  components  ua,  va,  and  wa  are  determined,  and 
subsequently  the  chaff  aerodynamic  forces  and  moments.  The  differen- 
tial equations  of  motion  comprise  only  body  velocities  and  accelerations 
resulting  from  these  aerodynamic  forces  and  moments. 

\ 

For  computational  purposes,  the  time  differential,  At,  appearing 
in  Equations  (D-22a,b,  c)  is  taken  as  the  integration  interval.  To  avoid 
an  indeterminancy  at  the  start  of  the  integration  process,  the  turbulence 
velocities  are  not  computed  until  after  three  Runge-Kutta  integration 
steps  (used  to  start  the  numerical  integration  of  the  equations  of  motion) 
are  completed.  Once  initiated,  the  time  dependency  of  the  turbulence  is 
forced  to  be  singular,  such  that  no  ambiguity  exists  if  a re-start  of  the 
main  predictor-corrector  integration  scheme  occurs. 
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The  following  NAMELIST  format  variables 
the  computer  data  inputs  for  j argument  values  (j 

*w  ■ WM  (J) 

3 Xu 

3h“  n W>MDZ  (J) 

Lz  = WLSCL  (J) 

0o  * WSIG  (j) 

which  are  inputs  at  the  altitude: 

h = TABZ  (J) 
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are  in  orporated  in 
* 10): 


APPENDIX  E 


A TV  OSPHERIC  TURBULENCE  DATA 


Generally  speaking,  turbulence  can  be  expected  in  a fluid  when 
there  is  shearing  flow  and  when  inertial  effects  are  much  larger  than 
viscous  effects.  In  the  atmosphere,  the  areas  of  principal  wind  shear 
are  near  the  ground,  in  jet  streams,  storm  cells,  and  clouds.  Secondary 
shear  layers  are  caused  and  influenced  by  geostrophic  winds  and  thermal- 
buoyancy  inequalities.  There  is  an  enormous  quantity  of  measured  data 
on  these  various  situations. 

Since  this  effort  is  concerned  with  "representative"  behavior, 
an  approximate  model  was  formulated.  Total  meteorological  detail  is 
not  within  the  scope  of  this  effort;  for  example,  non-isotropic  flow 
behavior,  such  as  jet  streams,  clouds,  etc.,  are  not  considered  (see 
Appendix  D,  Equation  (D-14)). 

Two  generalized  situations  are  discussed  in  this  appendix.  In 
Section  1,  single  valued  variables  define  the  free-atmosphere  turbulence, 
i.  e. , for  altitudes  above  the  earth  boundary  layer.  * Section  2 describes 
several  turbulence  parameters  as  a function  of  altitude  from  sea-level 
to  40,  000  feet. 

1.  Single  Valued  Turbulence 

The  relevant  data  to  the  digital  simulation  equations  are  the 
principal  scale  of  turbulence,  L,  and  the  standard  deviation  of  turbulent 
velocity,  oQ  . Reference  29  reviews  the  values  of  these  two  parameters 
as  used  in  various  theoretical  turbulence  ("frozen-flow")  programs. 
Reference  30  gives  a comprehensive  data  summary  of  medium- altitude 
turbulence  measurement  (MEDCAT).  The  principal  dimensional  values 
of  turbulence  obtained  from  these  and  other  references  are  presented 
in  Table  E-l.  The  power  spectral  functions  listed  in  this  table  are  de- 
scribed in  detail  in  Reference  30.  All  the  spectral  functions  in  the  table 
are  similar,  with  the  exception  of  that  from  Reference  31.  Reference 
31  incorporates  the  Von  Karman  spectrum,  generally  accepted  as  the 
most  accurate  distribution  function,  particularly  at  the  higher  wave 
numbers.  But  this  spectrum  is  a mathematically  irrational  function. 


* The  earth  boundary  layer  extends  to  2000  - 5000  feet  above  ground 
height. 
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Table  E-l. 


COMPILED  DATA  ON  MEDIUM  ALTITUDE 
TURBULENCE  PROPERTIES 


References 

- L - 
(Ft) 

°o 

Ft/Sec 

Power  Spectral 
Functions 

22 

6000 

6.9 

Von  Karman 

23 

5000 

1.0 

Dryden 

24 

1000-2000 

6.  0 

Dryden 

21 

1000-2000 

3.2 

Taylor-Bullen 

25 

1000 

1.0 

Press -Meadows 

20 


200-300 


3.2 


Dryden 


and  so  the  simpler  Dryden-type  spectrum  is  used  here  (the  actual 
numerical  difference  between  the  two  is  small).  The  Dryden  principal 
and  secondary  spectral  distribution  functions  are  respectively: 


**  (ft) 

*n 


2°o  LZ 
IT 


[_1  + (Lz  ft)2 

1 + 3 (Lz  ft)2 

(1  + (LZ  ft)2] 


(E-l) 


(E-2) 


The  spacial  correlation  function  corresponding  to  Equation  (E-l)  is 
Equation  (D-7)  in  Appendix  D. 

The  choice  of  the  scale  parameters  was  based  principally  on 
the  data  in  References  29  and  30.  Reference  31  relates  that,  "Gust 
Models  based  on  high-altitude  turbulence  are  in  general  found  to  be  too 
severe  and  unrealistic.  Based  upon  pilot  reactions  from  both  fixed-base 
and  flight- simulation  studies,  much  lower  values  of  gust  severity  and 
integral  scale  are  indicated  relative  to  those  often  quoted.  Values  of 
3 ft/sec.  for  rms  severity  ( aQ  ) instead  of  6 or  8 ft /sec.,  and 
L = 200  - 300  ft.  rather  than  2500  ft.  , were  found  to  be  much  more 
realistic.  " 


Reference  30  discusses  the  nonhomogeneous  character  of  the 
real  atmosphere  in  which  turbulence  is  noi  uniformly  distributed  but 
concentrated  in  pockets  (similar  to  CAT).  The  turbulent  region  is 
composed  of  two  sections  of  Gaussian  turbulence:  a high -amplitude 
section  which  consists  of  10  percent  of  the  region  with  an  rms  value  of 
8 ft/sec,  and  a low-amplitude  section  made  up  of  the  other  90  percent 
with  an  rms  value  of  2.  The  entire  region  has  a single  value  of 
0Q  i:  3.  16  ft/sec.  for  normal  fit,  which  is  quite  typical  of  real  turbu- 
lence. Figure  E-l  shows  an  averaged  model  proposed  for  this  situation 
in  Reference  30  , p.  139. 


Figure  E-2  shows  normalized  spectral  distributions,  (from 
Equation  (E-l))  for  three  scale  lengths  L,  and  compares  them  with 
normalized  data  from  Reference  30  , p.  94.  The  amplitude  of  the 
experimental  power  spectral  data  is  not  precisely  known  because  the 
spectrum  is  not  measured  to  all  wavelengths.  A common  method  of 
correlation  (Reference  30,  p.  217)  involves  laying  the  observed 
spectrum  over  the  theoretical  curves  and  obtaining  a match  near  the 
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Probability  of  Occurrence 


0 10  20 


Gust  Amplitude  - ft/sec 


Figure  E-l.  Theoretical  Mixture  of  Gaussian  Amplitudes 
of  Medium  Altitude  Atmospheric  Turbulence 


Normalized  Power  Spectrum,  $ - (cpf) 


/'  'A 
k\  / ' 


Flight  229,  Run  1 

Oy  ■ 3.3  ft/sec 
0y  * 4.4  ft/sec 
ow  » 3.1  ft/sec 


Dryden  longitudinal 
formulation 


K - Wave  Number  - cpf 


Figure  E-2.  Power  Spectral  Density  of  Turbulent  Spectrum 
Comparing  Dryden  Formulation  to  MEDCAT  Data 


"knee"  area.  Figure  E-2  indicates  i good  correlation  of  slopeo  on 
both  sides  of  the  "knee".  However,  note  that  this  case,  with  its  differ- 
ent values  for  ( oy,  oy»  and  o z)  is  not  truly  isotropic. 

For  the  most  generalized  situations  above  the  earth  boundary 
layer  the  following  values  are  proposed  for  chaff  dynamic  modeling: 

0Q  * 3.  Oft/sec. 

Lz  500  ft. 

Xy  £ 40  ft/sec. 

2.  Turbulence  as  a Function  of  Altitude 


The  integral  scale  Ljr  and  the  rms  turbulence,  oQ  , are  vari- 
ables with  altitude.  However,  above  the  earth  boundary  layer,  the 
Integral  scale  is  generally  assumed  constant.  The  rms  turbulence  can 
be  determined  from  the  expression: 


where  i is  the  turbulent  intensity  and  Xw  is  the  mean  wind  velocity. 
Reference  35  discusses  turbulent  intensity  in  the  free  atmosphere, 
recommending  a nearly  constant  value  above  the  ground  shear  layer  of; 

i = 0.  01  to  0.  03 

Figure  E-3,  obtained  from  Reference  35,  displays  the  uniform  turbulent 
intensity  above  the  ground  layer  (h  = 2000  to  5000  ft).  The  rms  turbu- 
lence intensity  is  generally  much  larger  within  the  boundary  layer  than 
within  the  free  atmosphere,  and  it  increases  with  ground  proximity. 
iBotropic  flow  does  not  exist  in  the  boundary  layer  where  the  rms  turbu- 
lence values  are  different  in  the  X,  Y,  and  Z directions.  Reference 
29,  p.  19  expresses  rms  turbulence  in  the  boundary  layer  in  terms  of 
shear  parameters.  All  three  components  are  assumed  to  be  identical 
with  the  value: 


oo  = 0.8  Xw/ln  (h/ZG)  (E-3) 

where  Zc  is  a height  scale  related  to  ground  protuberances  (typically 
1/30  the  characteristic  dimension  of  the  protuberance)  chosen  as  0.  116 
feet.  Wind  velocity,  Xw  , for  the  neutrally  stable  condition  (Reference 
36)  is  calculated  using  the  standard  1/7  power  law: 
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Figure  £-3 


The  Nominal  Intensity  of  Turbulent 
Fluctuations  as  a Function  of  Altitude 
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1/7 


(E-4) 


The  edge  of  the  boundary  layer  is  as  Burned  to  be  at  an  altitude  of 

= 40  ft/sec. 


hD  = 5000  ft,  with  a wind  velocity  of  X 


w. 


The  integral  length  scale  Lz  also  varies  within  the  earth 
boundary  layer.  Figure  E-4  compares  the  integral  scale  length  varia- 
tions with  altitude.  The  orders  of  magnitude  of  change  reflect  the 
radical  differences  between  the  ground  layer  and  the  free  atmosphere. 
References  37  and  38  express  the  length  scale  ao: 


LZ  = 0. 4 h 
LZ  = 2.  1 h' 


0.  73 


(h  < 465  ft)  (E-5) 

(465  ft  < h < 1 *00)  (E-6) 


Table  E-2  summarizes  the  turbulence  parameters  incorporated 
into  the  chaff  dynamic  model  for  preliminary  studies,  as  obtained  from 
the  previously  discussed  equations.  These  turbulence  values  are  only 
approximate  representations  • the  fundamental  variability  of  atmos- 
pheric flow  precludes  modeling  every  situation.  The  use  of  specific 
statistical  data  instead  of  this  generalized  model  is  recommended  when- 
ever experimental  and  theoretical  correlation  of  dynamic  behavior  or 
dispersion  is  being  performed. 
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Table  E-2.  TURBULENCE  PARAMETERS  WITH  NEUTRALLY 
STABLE  CONDITIONS  USED  IN  THE  COMPUTER 
NAMELIST  INPUT 


( TABZ) 
Ft 


5,000 


30.000 

40.000 


(WM) 

Ft/Sec 


(WSIG) 

Ft/Sec 


(WLSCL) 

Ft 


I 3h  | 
(WDMDZ) 
Ft/Sec/Ft 


. 0377 
.0084 


* a0  s 0.  03  Xw 

( ) parenthetical  quantities  are  NAMELIST  variables 
as  described  in  Appendix  D. 
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